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This  paper  develops  in  detail  the  theory  of  vorticity  generation  by  asymmetries  which  3re 
present  when  a  cylindrical  pulse  of  energy  is  deposited  in  a  gaseous  medium.  We  identify  three 
classes  of  asymmetry  which  may  occur  separately  or  in  combination,  depending  on  the  experiment 
being  analyzed.  We  also  introduce  a  convenient  representation  of  the  flow  field  in  terms  of  one  or 
more  vortex  filament  pairs  having  identical  strengths.  The  equations  which  give  the  vortex  filament 
strength  for  the  various  asymmetry  classes  differ  only  in  the  value  of  a  form  factor,  which  contains 
information  regarding  geometric  effects  and  hydrodynamic  interactions.  The  form  factor,  however.  . 

(Continued) 


00  ,  1473  EOiTiOn  Of  '  NOV  •>  I 


it  nktm  ri 


security  classification  of  this  page  Do  to  t*to*od) 


ABSTRACT  (Continued) 

depends  only  weakly  on  the  profile  of  the  density  distribution  of  the  medium  prior  to  energy 
deposition.  Finally  we  verify  and  calibrate  the  analytic  results  by  providing  data  from  accurate 
two-dimensional  simulations  of  several  cases.  The  analytic  and  numerical  calculations  agree 
closely  for  all  examples  considered. 
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VORTICITY  GENERATION  BY  ASYMMETRIC  ENERGY 
DEPOSITION  IN  A  GASEOUS  MEDIUM 


<c. 


I.  INTRODUCTION 

Recent  experiments  by  Greig  et  al.1-5  have  studied  hut  channels  produced  by  lasers  and  electric 
discharges  in  ambient  air.  Following  expansion  to  pressure  equilibrium,  the  hot  channels  cool  with  a 
rate  which  is  several  orders  of  magnitude  faster  than  that  attributable  to  classical  thermal  conduction. 
Since  any  reasonable  estimate  of  background  gas  velocities  falls  far  short  of  explaining  the  cooling  rate, 
we  require  a  mechanism  to  convert  curl-free  fluid  expansion,  which  does  not  mix,  into  persistent  vorti- 
city,  which  can.  The  importance  of  determining  the  mechanisms  and  scaling  laws  governing  the 
dynamics  and  cooling  of  such  hot  channels  derives  ftom  the  application  of  discharge  physics  to  other 
areas,  such  as  beam  physics  and  the  study  of  r.itrogen  fixation  in  the  atmosphere  by  lightning.'  In  the 
latter  case,  the  cooling  rate  of  the  hot  channel  gas  is  critical  to  predicting  the  global  production  of  nitro¬ 
gen  oxides. 

In  Section  II  of  this  paper  we  propose  a  mechanism  for  the  generation  of  persistent  Hows  which 
mix  cold,  ambient  gas  into  the  hot  gaseous  channels  produced  by  one  or  more  pulses  of  energy.  The 
mech  .nism  relies  on  deviations  from  cylindrical  symmetry  in  energy  deposition  by  a  given  pulse,  and 
we  identify  three  classes  of  asymmetry  which  appear  to  be  physically  significant.  Section  II  also  intro¬ 
duces  a  convenient  representation  of  the  vorticity  distribution  in  terms  of  one  or  more  vortex  filament 
pairs.  We  use  this  representation  in  Section  III  to  derive  a  formula  for  the  mixing  time  scale  rmiv 
which  characterizes  the  rate  at  which  cool  ambient  air  is  entrained  Into  the  hot  channel.  To  estimate 
rmix,  we  use  dimensional  analysis  to  derive  an  approximate  formula  for  ihe  expansion  induced  vorticity. 
Our  estimate  shows  that  the  proposed  mechanism  causes  mixing  of  cool  background  gas  with  the  heated 
channel  interior  on  time  scales  which  are  orders  of  magnitude  shorter  than  those  characterizing  molecu¬ 
lar  thermal  conduction.  Section  IV  provides  a  detailed  theoretical  analysis  of  the  three  classes  of  asym¬ 
metry.  The  resulting  formulae  for  the  vorticity  strength  are  considerably  more  accurate  than  that  of 
Section  III.  Section  V  presents  the  results  of  detailed  two-dimensional  simulations  of  several  sample 
cases.  These  computations  validate  the  qualitative  features  predicted  by  the  analytic  model  and  permit 
us  to  calibrate  the  analytic  formulae  accompanying  the  model.  We  have  used  our  results  to  analyze 
experimental  data  on  the  cooling  of  electric  discharge  channels  '  Our  comparisons  have  indicated  good 
agreement  between  predicted  and  observed  mixing  times. 
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II.  MECHANISMS  FOR  CONVECTIVE  COOLING 

The  symmetric  expansion  of  a  cylindrical  channel,  produced  by  a  pulse  of  energy  causes  no  con¬ 
vective  mixing  per  se.  but  asymmetries  in  energy  deposition  which  are  inherent  in  the  pulse  structure 
or  which  occur  relative  to  the  local  density  distribution  generate  long-lived  vortex  filaments.  The 
sheared  flow  attending  the  motion  and  interaction  of  these  filaments  mixes  cold  background  gas  into 
the  heated  channel  at  a  rate  which  depends  on  the  strength  of  the  induced  vorticity  and  thus  ultimately 
on  the  fluid-dynamic  asymmetries.  We  can  identify  three  generic  types  of  asymmetries:  ( I )  two- 
dimensional  asymmetries  from  pulse  displacement  off  the  axis  of  an  existing  hot  channel.  (2)  two- 
dimensional  distortions  of  the  pulse  envelope  from  a  circular  cross  section,  and  <3>  three-dimensional 
distortions  <e.g.,  curvature)  of  the  envelope,  such  as  characterize  a  lightning  or  spark  channel. 

As  the  channel  of  the  most  recent  pulse  expands  to  pressure  equilibrium,  any  deviations  from 
cylindrical  symmetry  will  lead  to  asymmetries  between  the  gradients  of  the  pressure  and  density  distri¬ 
butions.  The  equation  which  describes  the  resulting  vorticity  distribution  and  evolution  is 

df, 

f  V  v  -  f  •  Vv  +  (Vp  x  \7P)/p2  (1) 

where 


y  (2) 

is  the  vorticity,  _v  is  the  fluid  velocity,  p  is  the  density,  and  P  is  the  pressure.  All  of  the  variables  are 
functions  of  the  position^;  and  the  time  t.  Following  expansion  of  the  channel  to  achieve  pressure 
equilibrium,  significant  residual  vorticity  exists.  This  vorticity  is  responsible  for  mixing  ambient  gas 
with  the  hot  channel  gas.  Here  we  should  note  that  other  mechanisms  for  generating  residual  vorticity 
in  the  heated  channel  may  exist.  As  the  channel  expands  to  pressure  equilibrium,  for  example,  the 
Rayleigh-Taylor  instability  could  grow  significantly.  The  experimental  data  do  not  support  this  possibil¬ 
ity,  however.5  Another  possibility  is  the  rapid  movement  of  the  discharge  current  axis  which  could 
result  from  magnetic  forces  present  when  the  current  is  nonnegligible.  This  could,  in  turn,  displace  the 
surrounding  air  sufficiently  to  produce  some  long-term  mixing  motion.  We  do  not  attempt  to  estimate 
the  magnitude  of  such  phenomena  in  this  paper. 

Figures  1-3  depict  the  classes  of  asymmetry  schematically  and  will  form  the  basis  for  our  notation 
and  calculations,  as  discussed  in  subsequent  sections.  While  ail  three  types  will  often  occur  simultane¬ 
ously,  we  will  treat  them  separately  to  isolate  their  respective  characteristics.  We  have  defined  the  r- 
axis  to  be  collinear  with  the  channel  axis  (i.e.,  outward  from  the  page  to  the  reader)  for  cases  1  and  2. 
For  our  analysis  of  three-dimensional  distortions,  we  have  chosen  a  section  of  the  channel  with 
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moderate  sinusoidal  curvature  relative  to  the  length.  The  c-axis  is  collinear  with  the  center  of  the  mass 
of  the  section.  The  displacement  T0  represents  the  degree  of  asymmetry  which  is  present  in  a  given 
situation. 

Figure  1  shows  how  we  model  the  first  asymmetry  class,  which  is  relevant  to  sequences  of  approx¬ 
imately  collinear  pulses  and  to  pulses  with  interior  "hot  spots."  Our  model  assumes  a  preformed,  hot, 
low  density  channel  in  which  the  next  pulse  is  propagating  off  center.  In  Fig.  2  we  presen*  the  cross 
section  of  the  pulse  as  elliptical  where  the  ellipticity  is  defined  by  the  parameter  -Y0.  Residual  vorticity 
will  be  present  whenever  X0  is  nonzero.  Such  a  distortion  of  a  pulse  may  occur  in  several  ways:  ( 1 ) 
the  envelope  of  the  pulse  undergoes  smooth  spatial  oscillations;  (2)  the  pulse  is  deposited  in  a  region 
with  an  azimuthally  nonuniform  local  density  distribution;  or  (3)  the  source  of  the  pulse  produces  an 
azimuthally  nonuniform  (though  smooth)  distribution  of  energy.  Figure  3  shows  a  simple  model  of  the 
last  asymmetry  class,  which  includes  all  pulses  with  a  curved  axis.  Such  curved  discharges  occur  both 
in  laboratory  discharges  and  in  lightning.5  We  have  found  that  all  three  types  of  asymmetry  will  pro¬ 
duce  a  sizable  nonzero  value  for  which  will  cause  mixing  in  the  (r,  9)  plane.  In  the  case  of  three- 
dimensional  distortions,  will  also  be  nonnegligible;  this  results  in  mixing  in  the  (r,  :)  plane. 

III.  MIXING  TIME  SCALE 

A.  Vortex  filament  Representation 

We  may  represent  the  residual  vorticity  distribution  in  terms  of  one  or  more  vortex  filament  pairs 
of  strength  ±  ifla,  where  the  index  /  labels  the  asymmetry  class  and  a  denotes  z  or  9.  The  quantity  *,Q 
is  the  integral  of  the  vorticity  over  the  domain  containing  a  given  filament.  For  example,  in  the  case  of 
a  pulse  which  is  noncollinear  with  a  preformed  hot  channel  (asymmetry  type  1)  we  have  the  (x,  v)  flow 
pattern  depicted  in  Fig.  4.  These  flows  are  equivalent  to  those  of  a  single  vortex  filament  pair  having 
strengths  of  ±ku  and  respective  locations  (x,  ±.v).  The  strength  if,,  is  given  by  the  integral  of  the 
vorticity  over  the  upper  half  plane,  i.e., 

»f1;(r)  —  dy  f  ^  dx  f .  (x.  y,  r)  (3) 

and  the  coordinates  are 

x (t)  -  — - —  f  dy  f  dx  x  f.(x.  v.  r)  (4) 
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The  quantity  r  is  the  time  interval  over  which  vorticity  generation  is  completed  and  is  thus  approxi¬ 
mately  me  time  required  for  the  hot  channel  gas  to  expand  to  pressure  equilibrium.  Defining  t  -  0  to 
be  the  time  at  which  the  discharge  is  initiated,  we  may  integrate  Eq.  (1)  over  the  interval  <0,  r)  to 
obtain  the  vorticity  distribution  €,  (x.  y,  r). 


Figure  2  shows  that  the  elliptical  channel  will  have  two  associated  vortex  filament  pairs.  Here  the 
domain  of  integration  in  Eqs.  ( 3 ) - ( 5)  will  be  the  quarter  plane  and  we  will  use  cylindrical  coordinates. 
For  the  curved  channel  section  in  Fig.  3,  we  will  again  use  cylindrical  coordinates.  The  integral  for  Hiz 
will  cover  the  upper  half  plane  as  in  Eqs.  (3)-(5);  however,  unlike  class  one,  *i:  will  vary  along  the 
channel  axis.  To  obtain  we  will  integrate  over  half  a  wavelength  in  our  simple  model  of  a 
sinusoidally  curved  channel.  Because  will,  in  general,  vary  with  9 ,  we  may  then  perform  an  azimu¬ 
thal  integral  (0«[O.  ir])  to  obtain  a  total  longitudinal  mixing  strength. 


B.  Derivation  of  Mixing  Time  Scale 

In  order  to  derive  an  equation  for  the  time  scale  of  the  mixing  of  cool  ambient  gas  with  the  hot 
channel  gas,  we  must  first  review  the  experimental  observations.1-5  Soon  after  the  deposition  of  energy 
(=300  -  600  J/m)  by  an  electric  discharge  or  laser  pulse,  the  hot  channel  expands,  producing  a  shock 
wave.  Within  100  ^s,  the  shock  has  decoupled  from  the  channel,  which  has  reached  pressure  equili¬ 
brium  and  ceased  to  expand.  At  that  point,  the  channel  boundaries  are  smooth.  However  after  another 
100  ms,  the  edges  of  the  channel  have  become  distorted,  and  the  channel  has  started  to  expand  and 
coo!  'hrough  entrainment  of  the  surrounding  air.  By  =  1  ms.  small  scale  (turbulent)  structure  is  evi¬ 
dent  and  the  channel  has  grown  considerably.  We  attribute  the  initial  disruption  of  the  channel  and  the 
entrainment  of  ambient  air  to  the  convective  motion  generated  by  the  mechanism  of  Eq.  (1).  To  deter¬ 
mine  the  mixing  time  scale,  therefore,  we  will  treat  this  large  scale  structure  separately  from  the  smaller 
scale  motion  observed  as  the  channel  cools.  An  understanding  of  the  development  of  the  latter  flow 
structure  would,  of  course,  be  necessary  for  an  accurate  description  of  the  dynamics  of  the  fluid  inside 
the  channel. 

For  simplicity  we  assume  that  the  asymmetry  induced  flow  field  may  be  represented  approximately 
as  two  compact  vortex  filaments  of  strength  ±i?  at  ±y  ,  as  in  Fig.  4;  thus  we  treat  each  vortex  filament 
pair  separately.  This  is  in  fact  the  situation  for  f  .  in  the  case  of  off  center  beam  propagation.  The 
azimuthal  velocity  induced  by  each  filament  decays  as  1/r  away  from  the  vortex  center.  The  velocity  of 
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the  fluid  along  the  symmetry  plane  is  the  sum  of  flows  induced  by  each  vortex  separately.  Figure  4  is  a 
schematic  of  the  situation  with  a  sketch  of  the  variation  of  the  fluid  flow  velocity  along  the  x-axis.  This 
velocity  is  given  by 


vt(jc,  0) 


KV 

rrU’  +  v2) 


The  filaments  are  migrating  in  the  same  direction  at  a  slower  velocity, 


(6) 


vv  -  17) 

which  we  will  ignore  in  the  following  integral  estimates.  This  overall  migration  velocity  can  be  impor¬ 
tant  when  a  series  of  pulses  is  being  considered.  In  that  case  a  quasi-steady  state  develops  in  which  the 
systematic  migration  of  the  integrated  vorticitv  entrains  cold  fluid  stochastically  at  the  edges  of  the  hot 
channel. 


To  estimate  the  mixing  time,  we  use  Eq.  (6)  to  calculate  the  time  required  for  a  fluid  element 
starting  at  x  -  -S0  to  reach  x  -  50,  ignoring  the  effects  of  small  scale  turbulence.  The  equation  of 
motion  for  this  element  is 


dx 

dr 
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which  can  be  integrated  to  give 


The  quadrature  can  be  performed  and  gives  the  following  mixing  time  estimate: 
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where  the  vortex  strength  ic  is  given  by  Eq.  (3)  and  the  "average"  displacement  of  this  vorticity  above 
the  y  -  0  plane  is  given  formally  by  Eq.  (5).  Recently  we  have  found  that  rmix  has  a  simpler  and 
perhaps  more  realistic  interpretation.  If  we  define  V0  to  be  the  volume  of  the  channel  just  after  pres¬ 
sure  equilibrium,  rmix  is  approximately  the  time  interval  required  for  a  volume  1'0  of  cool  ambient  air 
to  mix  with  the  hot  channel  gas.  For  an  ideal  gas,  therefore,  rmix  is  approximately  the  time  in  which 
the  channel  volume  doubles.6 


We  note  in  Fig.  2  that,  for  the  case  of  the  elliptical  channel  (asymmetry  class  2),  two  vortex 
filament  pairs  are  generated  and  the  flow  is  inward  along  the  semimajor  axes,  at,  and  outward  along  the 


semiminor  axes,  bt.  A  representative  rmi<  (and  volume  doubling  time)  can  then  be  taken  as  the  sum 
of  the  time  intervals  required  for  a  fluid  element  to  traverse  half  the  ellipse  along  a,  and  be  respec¬ 
tively.  For  these  times  and  each  vortex  filament  pair,  the  derivation  of  Eq.  (10a)  will  carry  over 
directly.  For  the  inward  flow  along  the  semimajor  axis,  a, ,  we  integrate  from  0  to  ae  for  the  fluid  parti¬ 
cle  path  rather  than  -S0  to  S0.  In  Fig.  2,  the  vortex  filament  pair  most  influencing  this  motion  is  that 
appearing  on  the  right  side  of  the  minor  axis.  The  displacement  of  these  filaments  is  2 y  =  2 be.  This 
gives  us 
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Similarly,  for  a  fluid  element  traveling  along  the  semiminor  axis  from  y  -  0  to  be,  the  vortex  filaments 
most  influencing  the  motion  are  the  two  in  the  upper  half  plane  and  are  therefore  displaced  by 
2.t  =  2a, ;  so  we  have 


1  mix<b) 


nbe 
Ik2.  la. 


(10c) 


Our  representative  time  scale  would  then  be  the  sum  of  Eqs.  (10b)  and  (10c), 


rrmx(2(  **  '  mtxU)  '  rnix'bl 


( 10d) 


The  case  of  three-dimensional  distortions  is  more  complex  than  the  others,  since  two  components 
<kj.  and  k *)  of  the  residual  vorticity  are  nonnegligible  and  are  functions  of  c  and  9,  respectively.  From 
Fig.  3,  we  see  that  for  a  given  value  of  r.  it,,  is  computed  by  assuming  that  the  density  distribution  is 
offset  from  the  pressure  distribution,  as  for  asymmetry  class  1.  The  equation  for  is  therefore 

identical  to  Eq.  (10a)  with  k  as  a  function  of  r.  We  will  find  that  the  dependence  on  :  leads  to  localiza¬ 
tion  of  the  vorticity  at  r  •  0,  \/2,  and  in  contrast  to  the  case  of  off  center  beam  propagation.  Simi¬ 
larly  we  will  find  that  !k*I  is  largest  when  9  -  0°  or  180°,  so  that  we  can  compute  as  the  time 

for  a  fluid  element  to  travel  half  a  wavelength  (for  a  sinusoidally  curved  channel)  under  the  influence 
of  vortex  filaments  of  strength  !k*(0°)I  -  |k*(180°)I.  With  a  separation  2.v  =2 R{.  where  R\  is  the 
radius  of  the  channel  cross  section  at  pressure  equilibrium  we  obtain 


7 r\ 

2/?,  Ik*  <0)1 


(10e) 
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While  Eqs.  (10a)  and  (lOe)  will  give  us  an  estimate  of  the  mixing  rate  for  asymmetry  class  three, 
we  can  define  a  channel  doubling  time  only  by  computing  the  flux  of  ambient  air  into  the  channel 
under  the  influence  of  the  total  velocity  field  defined  by  the  strengths  ic3.(r)  and  #< ( « > . 

C.  Estimates  of  rmU  from  Dimensional  Analysts 

We  may  now  use  Eq.  (10a)  to  estimate  rm,x  for  the  first  asymmetry  class,  in  which  a  pulse  is  not 
collinear  with  a  preformed  hot  channel.  To  do  so.  we  require  values  for  ic  and  y.  The  remainder  of 
this  paper  deals  with  the  analytic  calculation  and  numerical  calibration  of  expressions  for  ic.  We  will 
find  a  reasonable  value  for  F  from  the  numerical  simulation  of  sample  problems.  A  simple  dimensional 
analysis  allows  us  to  make  at  least  crude  estimates  rather  directly  for  the  example  of  two  successive, 
noncollinear  pulses.  A  number  of  size  scales  enter  the  problem:  S0.  the  radius  of  the  channel  cr  ^ 
by  the  first  pulse;  the  characteristic  scale  lengths  for  the  pressure  and  density  gradients;  and  the  *1 

and  final  radii  R»  and  R)  of  the  second,  displaced  channel  being  formed.  Fortunately  most  of  e 

scales  are  either  unimportant  or  expressible  in  terms  of  S0.  If  the  first  and  second  pulse  have  s 
initial  overpressures,  R  |  and  So  will  be  roughly  equal  and  R$  will  be  a  modest  fraction  of  R\.  .a 

numerical  simulations  (Section  V)  support  the  approximate  equality  of  S0  and  Rt. 

W'e  expect  other  simplifications.  W  hether  a  shock  expansion  or  an  adiabatic  expansion  is  being 
driven  by  energy  deposition,  the  larger  the  pressure  gradient  the  smaller  the  time  over  which  it  acts. 
Thus  the  integrated  vorticity  is  relatively  insensitive  to  the  shock  thickness.  Similarly,  the  density  gra¬ 
dient  is  integrated  over  space  so  the  inner  and  outer  densities  fp0  and  p*)  enter,  but  the  scale  length 
of  the  transition  region  can  be  neglected.  W'e  expect  and  will  in  fact  show  that  the  maximum  vorticity 
generation  occurs  when  ,V0  =  S0.  Thus  there  are  no  small  parameters  arising  as  ratios  of  characteristic 
lengths  unless  we  consider  exceltionally  tightly  focused  pulses,  in  which  case  /?p  «  R\  =  S0,  or  only- 
slight  departures  from  the  symmetric  superposition  of  the  pulses,  in  which  case  Xn  «  S0. 

The  integrated  vorticity  k  has  units  of  cnr/sec.  a  length  times  a  velocity  or  a  characteristic  time 
multiplied  by  a  velocity  squared.  The  characteristic  velocity  will  be  the  expansion  velocity.  When  the 
energy  deposition  is  fast  and  the  pulses  are  strong,  this  velocity  is  a  characteristic  sound  speed  cs.  The 
characteristic  time  is  the  expansion  duration  (if  the  energy  deposition  is  slow)  or  the  sonic  transit  time 
S,y'cf.  Thus  we  expect 

k,.  =  c-j  S0  F i-  (11) 

where  F is  a  dimensionless  form  factor  containing  geometric  effects,  detailed  hydrodynamic  interac¬ 
tions.  and  information  about  the  channel  and  beam  profiles.  Cancellation  effects,  leav  ing  v  ■  V  v  as  the 
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dominant  contribution  to  the  vorticity  source  team,  generally  reduce  F somewhat  below  unity  W e 
will  demonstrate  this  in  Section  IV 


If  we  assume  that  y  =  0.85,).  Eq.  (10)  for  the  mixing  time  becomes 

rmix=  7.9S,y.  M2) 

For  a  strong  pulse  which  heats  the  gas  appreciably  on  passage,  we  shall  assume  F]:  =  0.5.  Taking 
5  x  104  cm/s  as  a  generic  value  for  cs  and  choosing  So  =  I  cm  gives  rmix  =  300  /as.  This  time  scale  is 
at  least  three  orders  of  magnitude  shorter  than  that  of  classicial  thermal  conduction.' 

If  we  consider  only  the  large  scale  convective  flow-s.  this  300  m s  "mixing  time"  from  the  dimen¬ 
sional  analysis  estimate  of  Eq.  (10a)  is  the  time  for  cold  material  from  one  side  of  a  hot  channel  to 
cross  the  channel  and  reach  the  other  side.  Because  of  the  presence  of  small  scale  turbulence  and  other 
perturbations  which  can  affect  the  interaction  of  the  vortex  filaments,  the  fluid  particles  will  be  likely  to 
follow  a  more  random  path  inside  the  channel.  In  addition,  true  mixing  of  hot  and  cold  fluid  will  prob¬ 
ably  take  somewhat  longer  because  one  or  two  rotations  of  the  vortices  will  be  required  to  entrain  and 
smear  in  an  appreciable  amount  of  the  cooler  fluid.  This  "mixing1  time,  however  best  defined,  is  clearly 
important  when  comparable  to  or  less  than  the  time  interval  between  pulses. 

In  the  next  section,  we  perform  an  integration  of  the  vorticity  source  term  in  Eq.  M)  to  improve 
on  Eqs.  (11)  and  (12),  our  dimensional  analysis  estimates.  Using  a  few  reasonable  assumptions,  vve 
obtain  a  quadrature  for  the  integrated  vorticity  k,q,  which  will  display  the  various  nondimensional 
dependences  of  the  form  factor  Fa. 

IV.  APPROXIMATE  VORTICITY  INTEGRALS 

Equation  (3)  defines  the  residual  vortex  strength  following  expansion  to  pressure  equilibrium 
lr  >  r)  and  will  be  evaluated  analytically  to  replace  Eq.  Ml)  in  the  mixing  time  estimates.  Eqs.  (10). 
First,  however,  we  must  integrate  Eq.  (1)  over  the  time  interval  (0.  r ) .  We  begin  by  stating  several 
assumptions  which  will  permit  us  to  perform  these  integrals. 

We  model  the  deposition  of  energy  as  instantaneous  and  use  pressure  pulses  of  finite  size  with  in 
appropriate  radial  profile  and  a  total  energy  which  is  equivalent  to  that  of  the  laser  pulses  or  electric- 
discharges.  This  is  a  reasonable  assumption  for  the  experiments  of  Greig  et  al.,1'  in  which  the 
discharges  were  much  shorter  in  duration  than  the  expansion  times  of  the  hot  channels.  Our  model 
deposits  the  energy  as  internal  energy  only,  and  we.  therefore,  do  not  consider  relaxation  of  excited 
molecular  and  atomic  states.  The  presence  of  long-lived  states  will  slow  expansion  of  the  channel; 
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however,  our  analytic  model  becomes  more  accurate  in  this  case,  as  discussed  below.  The  variable 
Ri"  will  denote  the  radius  of  an  expanding  cylinder  of  hot  gas  which  is  produced  by  the  latest  pulse. 
Note  that,  according  to  Section  III. B.  R(0)  -  Rt}  and  R<t)  *■  /?,.  To  represent  the  flow  field  of  the 
expanding  channel,  we  assume  that  the  flow  in  the  outside  region  ( r  <  Rit))  behaves  incompressibly 
and  the  inside,  beam-heated  region  expands  uniformly  Thus 


v,(r,  /) 


rUU)!R  it) 
i'U)  R  U)/ r 


if  r  ^  Rd) 
if  r  >  R  ([I 


(13) 


specifies  the  flow  everywhere  as  a  function  of  the  heated  region  radius  R  <r)  and  velocity  Lit)  =  Rit). 
This  flow  has  a  uniform  but  time  varying  divergence  inside  the  heated  region  and  zero  divergence  out¬ 
side.  In  reality  the  fluid  inside  the  just-heated  channel  will  give  up  energy  to  the  cold  surrounding  fluid 
via  shocks,  and  a  fraction  of  the  pulse-deposited  energy  will  even  escape  to  infinity  as  an  acoustic  wave. 
The  smooth  shape  of  the  expanded  channels  observed  from  ID  hole-boring  calculations'  and  the  com¬ 
putational  and  experimental5  result  that  most  of  the  deposited  energy  stays  close  to  the  original  pulse 
deposition  region  support  the  approximations  implied  by  using  the  flow  field.  Eq.  (13).  We  note  that 
Eq.  1 13)  becomes  a  better  representation  as  the  rate  of  energy  deposition  decreases. 


As  we  will  soon  show  more  explicitly,  the  important  feature  of  vorticity  generation  is  the  radial 
distance  which  each  fluid  element  moves.  Pressure  gradients  arising  from  accelerations  of  a  fluid  ele¬ 
ment  do  not  really  contribute  to  vorticity  in  the  present  context.  Because  the  fluid  elements  begin  and 
end  their  expansion-induced  displacement  at  rest,  the  average  acceleration  is  zero.  The_v  •  Vy  vorticity 
source  term  has  the  same  sign  throughout  the  expansion  and  consequently  contributes  more  strongly  to 
the  integrated  vorticity. 


Because  our  velocity  field  during  expansion  is  radial  and  varies  only  with  r.  the  coupling  term 
£  Vv  in  Eq.  (1)  is  initially  negligible.  Further,  the  maximum  vortex  induced  flow  speed  which  resides 
in  the  system  is  much  smaller  than  the  maximum  expansion  speed.  This  term,  therefore  remains  small 
relative  to  the  other  terms  throughout  the  expansion  and  may  be  neglected.  To  evaluate  the  pressure 
gradient  in  the  source  term  we  may  use  the  flow  field,  Eq.  (13).  in  the  equation  of  motion  to  obtain 


dt 


r-7- ( U/R)  +  r  U:/R2 
Jr 


-  4<U?  )  -  L  2/?2/r3 
r  dt 


i f  r  ^  Rit) 
if  r  >  R  it) 


(14) 
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Notice  that  we  assume  that  the  acceleration  in  the  radial  direction  has  values  strictly  appropriate 
only  when  feedback  of  the  asymmetric  density  gradients  on  the  driving  expansion  flow  are  small.  The 
model  appears  to  work  quite  well  for  large  density  variations  as  well,  a  result  that  is  understandable  in 
hindsight.  The  maximum  vortex-induced  flow  speed  wh'"h  resides  in  the  system  is  much  smaller  than 
the  maximum  expansion  speed.  Further,  the  vorticity  is  generated  essentially  instantaneously  relative 
to  the  mixing  timescale  rm„.  Thus  the  generation  term  can  be  calculated  assuming  that  the  asymmetric 
density  gradients  do  not  change  the  expansion-driven  pressure  gradients  and  that  the  vorticity  which 
develops  does  not  affect  the  density  gradients  during  the  relatively  brief  expansion. 

Since  the  pressure  gradients  are  assumed  to  be  radial,  the  radial  vorticity  component  will  be  negli¬ 
gible  and  we  need  to  evaluate  only  the  azimuthal  and  axial  gradients  of  the  density.  Our  representation 
of  the  density  profile  is 

P  (s.  r>  -  p„  exp  [-In  (p./po^  S  <£.  SG.  /)]  (15) 

where  .?<0,  S0,  0)  -  1.  g<°°.  So.  0)  -  0.  s  is  a  displacement  vector:  and  S0  is  a  characteristic  scale 
length  defined  in  Figs.  1-3.  Normally  g(j,  S0,  r)  is  a  function  of  the  ratio  _s/S0.  In  the  subsections 
below,  we  will  evaluate  kjo  for  the  three  general  asymmetry  classes.  For  brevity  we  will  discuss  the  off 
center  beam  propagation  case  in  detail  and  shorten  the  presentation  for  the  other  two  cases. 


A.  Off  Center  Beam  Propagation 


As  indicated  by  Fig.  1.  the  case  of  off  center  beam  propagation  is  two-dimensional,  and  the  dom¬ 
inant  contribution  is  the  passage  of  the  shock  across  the  preformed  hot  channel,  which  is  assumed  to  be 
cylindrical.  Thus  we  may  use  the  flow  field  at  and  outside  of  the  boundary  of  the  expanding  channel 
(r  >  R  (t)).  Because  the  problem  is  two-dimensional,  v.  p,  and  P  do  not  vary  with  r,  and  for  this  rea¬ 
son.  only  (.  is  nonzero.  Noting  that  the  flow  field  in  the  outer  region  is  incompressible  and  using  Eqs. 
<1>.  (14),  and  (15),  we  have  for 


dj. 

dt 


=  0 


In 

Poo 

dg(s) 

I1™1’  '  ‘'uJ 

Po 

SO 

r*  r2  dt 

if  r  <  R(t ) 
if  r  >  R(t) 


(16) 


Because  we  assume  that  the  vorticity  does  not  affect  the  density  gradients  during  the  relatively 
brief  expansion.  g(s)  does  not  vary  with  time  in  our  model,  and.  therefore,  describes  the  channel  prior 
to  deposition  of  the  latest  pulse.  Notice  that  we  have  suppressed  J0  as  an  argument  o f  .?(_£>■  W'e  "ill 
take  5,,  into  account  when  evaluating  form  factors  later  in  this  section.  We  may  now  calculate  icL.  by 


using  Eq.  <16)  in  Eq.  (3).  Since  we  are  currently  using  Eulerian  variables,  we  may  reverse  the  order  of 
the  time  ana  spatial  integra'.icns  to  obtain 


K1:(r)  -  In  <P«/Po>X,  dl  $;  d9  X,,,  dr  r 


Bgls) 


l  UR)2  _  _!_  d(UR) 

r*  dt 


We  now  change  to  the  Lagrangian  variables  <r0,  9I})  where 
r  —  r<t)  —  V  '’o  +  —  Rtf 


9  -  9  a. 


This  gives  us 


r  dr  =  r0  dr0, 


and  rh)  is  the  instantaneous  radial  position  of  a  fluid  element  initially  at  r < 0)  =  r0.  In  terms  of  rn, 
the  radial  density  variable  s  becomes 

s  -  (,Y0:  +  r2  -  2X0r0  cos  »)'  2  (20) 

and  the  integral  in  Eq.  (17)  can  be  performed  without  approximation.  This  gives  us 


r  d9 

Jo  Art 


Qg  giro  +  X0)  ~  glX o  —  r0),  r0  <  .Yq 
d9  ”  ^(r0  +  T0)  -  ?(r0  -  T0).  r0  >  X0 


and  Eq.  ( 1 7)  becomes 


Kl: 

1  \  1 
(r)  =  In  — 

P  0 

r 

Jr  a 

drarQ 

r2U) 

glr0 

+  X0) 

IUR): 

r2(t) 

1 

*1 

r* 

a  dr ,-jr a 

glX0 

~  r0) 

(UR)2 

d(UR) 

Jr 

r2(t) 

r2(t) 

dt 

r°° 

dr  ar  a 

(UR)2 

d(UR) 

r2(t)  gir°- 

T0) 

r2(r) 

dt 
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Notice  that  the  expression  in  Eq.  ( 22)  is  long  because  of  the  requirement  that  the  radial  variable  s 
be  positive  or  zero  Equation  (21)  is  necessary,  for  instance,  when  gis.  Sn)  is  a  square  well 

1 1 ,  s  ^  $0 

*<£■  Sol-jo.  s  >  S0  '  1231 

which  is  not  defined  for  negative  s.  However,  if  £<j.  S0)  is  an  even  function  of  s,  such  as  a  Gaussian, 
then  Eq.  (21)  represents  an  irrelevant  distinction  and  Eq.  (22)  shortens  considerably. 

The  integral  in  Eq.  (22)  is  difficult  to  perform  even  if  R  (r)  for  our  particular  case  is  known,  and 
further  approximations  must  be  made  to  get  a  usable  analytic  result.  The  approximation  most  useful  is 
to  evaluate  the  integrand  at  a  specific  time.  If  we  assume  that  the  initial  and  final  states  of  the  expan¬ 
sion  are  at  rest,  we  have  £/( 0)  -  U(r)  -  0.  The  expansion  flux  R(i)UU)  will  then  peak  at  a  time 
0<  r„  <  r.  We  will  evaluate  the  integrand  in  terms  of  rm.  If  the  function  R  <t)UU)  is  approximately 
symmetric  about  r  —  r„,  a  reasonable  approximation  for  the  term  L'2U)R2U)  in  Eq.  (22)  would  be  half 
the  maximum  value,  i.e.. 


U2it)R2U)  =  U2Um)R2UJi  2  =  U2R2! 2.  (24) 

We  also  replace  rtr)  by  rm  »  V rj  +  R2  -  R}t .  Because  rm  does  not  vary  with  time  and  because 

the  system  is  at  rest  at  r  =*  0  and  r  «  r.  the  term  proportional  to  integrates  to  zero.  This 

at 

term,  therefore,  corresponds  to  a  transient  in  the  vorticity  during  the  expansion  to  pressure  equili¬ 
brium.  As  the  time  approaches  t  -  r,  we  expect  the  vorticity  strength  to  approach  a  value  which  will 
be  approximately  constant  for  r  >  r.  W'ith  the  above  simplifications,  we  may  perform  the  time 
integration  of  Eq.  (22)  to  obtain 


k | .  (r )  ~ 


2 


In 


Poc 

Po 


r-  .  R^g{ro  +  x0) 
I**  r°'°W  +  Ri-  Ro) 


(25) 


~Sg!dr^ 


RmfflX  0  ~  '0> 
(r0:  +R2-  Rj  )2 


Jv  dr*rt) ' 


-  X0) 


W  +R2-  R2)2 


for  the  residual  vorticity.  We  may  approximate  the  expansion  time  r  by 


r  =  2(*,  -  R{))/Um 


(26) 


where,  as  before,  Ro  -  /?  ( 0)  and  R\  -  R(r).  In  Eq.  (26).  we  use  Vm/2  to  estimate  the  average  velo¬ 
city  of  expansion. 


Next  we  write  the  integral  in  terms  of  three  nc  ndimensional  parameters 


a  =  R„  <  1, 
b  =  Rm/Sc  ^  1, 
r  =  Xo/S  =  1. 


Letting  the  integration  variable  be  17  =  r/Rm  yields 


k  |_.  (r) 


d  v  y 


g  I50[i7^  +  c]) 
<1  -+-  —  a2)2 


g(S0[c  -  T)6l> 

(  1  -r  T)'  —  fl')’ 


gtSphb  -  c ] )  j 
(1  +TT-  a:):) 


127) 


(28) 


If  the  remaining  composite  integral  is  identified  as  a  three  parameter  form  factor  ft. (a,  b,  c)  with 
values  of  order  unity,  the  very  useful  improvement  on  the  dimensional  approximation  (Eq.  (ID)  is 
obtained: 

?|.  <t)  ~  < U2r/2 )  In  (pJpp)J\Aa.  b.  c> 

=  Um(R\  -  R n)  In  (pJpp)f\Aa,  b,  c).  (29) 

Here  the  form  factor  /,.  is  generally  less  than  1/2.  The  sign  of  f\  indicates  the  direction  of  flow  in  the 
(.v.  y)  plane— counterclockwise  or  clockwise. 


We  will  now  include  the  explicit  dependence  of  the  function  g(s.  50)  upon  S0;  usually  this  occurs 
as  the  ratio  s/S0.  In  this  case,  the  multiplicative  S0  in  Eq.  (28)  will  cancel  with  that  in  the  denomina¬ 
tor.  For  this  reason,  the  form  factor  j\.  in  Eq.  (29)  depends  only  on  the  ratios  a,  b,  and  c.  and  not 
explicitly  on  S0.  and  we  have 


f,.(a,  b,  c)  -  f  dt\r\ 


g  (r}b  +  c) 

(1  +  n2-  a2)2 


n 


g(c  -  rib) 

( 1  4-  r)2  —  a2)2 


g(yb  -  c) 

U  +  t,-  -  a2)2 ' 


(30) 


;  The  general  integral  form  factor,  Eq.  (30),  can  be  evaluated  numerically  for  any  reasonable  profile 

g(s/50).  Figure  5  shows  plots  of  J\.  versus  c,  the  measure  of  channel  separation,  for  several  values  of 

1 
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a  at  b  *»  0.7  and  for  several  values  of  b  at  a  —  0.6.  These  values  are  close  to  the  actual  initial  condi¬ 
tions  for  the  detailed  calculation  of  the  next  section.  In  this  figure  the  super-Gaussian  density  profile 
was  used. 


gsc(s/S0)  -  exp  <-  s2/S02) 


(31) 


and  the  integrations  were  performed  using  a  numerical  quadrature  algorithm.  We  have  also  evaluated 
/,.  for  a  Bennett  profile. 


Sfl(s/S0) 


Poo  ~  Pi) 

Pood  +  s2/sn2>2 


/  ln(p„/p0); 


(32) 


the  results  appear  in  Fig.  6.  In  Eq.  (32),  the  "Bennett  radius"  is  equal  to  S0.  In  the  case  of  a  square 
well  density  discontinuity  at  5  -  S0,  Eq.  (23).  the  integral  can  be  performed  analytically  to  give 


fit  (a.  b,  c ) 


J_  _ tr _ 

2  ( 1  +  c  )2  +  b2(  1  —  a2) 


_ b^ _ 

( 1  —  c)2  +  b2(  1  —  a2) 


(33) 


for  c  >  1  +  ab.  This  is  the  case  when  the  second  pulse  is  wholly  outside  the  channel  formed  by  the 
first  pulse.  As  c  approaches  °°,  / f?  decreases  as  —  1  /c3;  so  pulse  channels  separated  by  more  than  3 
radii  interact  only  weakly. 


When  the  channels  are  closer,  we  have  c  <  1  +  ab  .  and  the  channels  overlap  somewhat.  As 
long  as  c  >  1  -ab  there  is  still  part  of  the  second  channel  of  initial  radius  R0  outside  the  cylinder 
s  -  S0.  For  this  region  of  values,  1  -  ab  <  1  +  ab,  we  have 


/f7  ( a .  b,  c) 


2  ( 1  +  c);  +  b2(l  -  a2) 


(34) 


Equation  (34)  agrees  with  Eq.  (33)  at  the  interfacial  separation  c  —  1  +  ab. 


There  is  a  third  region  where  the  initial  radius  of  the  second  channel  lies  wholly  within  the  first, 
i.e.,  1  -  ab.  Then  <a,b,c)  satisfies  Eq.  (33)  again.  In  the  intermediate  region 

1  -  ab  <  c<  1  +  ab,  /j*  takes  on  its  lowest  negative  value  at  the  exterior  touching  point  c  -  1  +  ab. 
For  larger  values  of  c  the  magnitude  I  /j*  l  decreases  and  for  smaller  values  I  /]“  I  decreases  monotoni- 
cally  to  zero  at  c  -  0.  This  lowest  negative  value  is 


JT: 


^max 


-  2d  4-  ab) 
(4  +  4  ab  +  b2) 


(35) 


When  ab  -  b2  -  1/2,  i/f7  !max  -  6/13,  which  is  close  to  the  values  in  Figs.  5  and  6  for  /f?  and  /f.. 
respectively. 
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Also  of  interest  is  the  slope  at  c  *■  0,  where  the  pulses  are  concentric  and  no  vorticity  genera¬ 
tion  is  expected. 

-  -2b*  ,u» 

■o  (1  -  b2(\  -  a2)]2' 

which  is  nonzero  and  large.  Thus  even  modest  nonconcentricity  leads  to  appreciable  vorticity  and  mix¬ 
ing.  Similar  behavior  is  observed  in  the  two  smooth  profiles  considered  as  well.  Figure  7  compares 
ff:>  /f?.  and  /}*  for  a  —  b  —  l/%/2  to  demonstrate  the  similarity  of  the  vorticity  generation  form  fac¬ 
tors  regardless  of  the  profile.  This  relative  profile  insensitivity  arises  because  is  an  integral  quantity, 
the  total  upper  half  plane  circulation. 


a/K 

dc 


Returning  to  our  integrated  vorticity  estimate  of  Eq.  (29)  with  a  maximum  form  factor  of  6/13, 
and  Um  =  cs  we  have 


^R\  *0>  It! 


Po 


(37) 


Using  Eq.  (8)  and  assuming y  =  -  R0  gives  an  estimate  for  the  maximum  velocity  on  the  .x-axis: 


IvJmax  iJfl.  *n 


Po 


C'. 


(38) 


Under  optimum  conditions  the  velocity  between  the  vortex  filaments  approaches  a  quarter  of  the  sound 
speed  in  the  surrounding  fluid.  When  the  expansion  is  subsonic  throughout  (where  this  analysis  ought 
to  be  most  accurate).  cs  in  Eq.  (38)  should  be  replaced  by  the  appropriately  averaged  expansion  velocity 

Um- 


Using  Eq.  (10a)  for  the  mixing  time  withy  =  0.8S0  =  /?,  -  R0  and  it,.  Eq.  (37)  above  gives 

6.87tS0 

' -  (39) 

‘  IM 

With  c,  =  5  x  104  cm/s.  In  (pjp^)  =  2.5,  and  S0  “  1.0  cm.  the  fastest  mixing  time  is  about  170  ps. 
This  estimated  fastest  mixing  time  is  within  a  factor  of  two  of  the  dimensional  estimate.  In  any  particu¬ 
lar  system  the  integral  of  Eq.  (30)  can  be  performed  numerically  and  the  resulting  form  factor  substi¬ 
tuted  into  Eq.  (29)  for  the  integrated  vorticity. 

The  sign  of  is  negative  for  the  configuration  of  Fig.  1  where  the  second  pulse  is  centered  to 
the  left  of  the  original  channel.  This  means  that  the  jet  of  colder  gas  across  the  original  channel  starts 


In 
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on  the  side  opposite  to  the  displaced  second  pulse  and  rushes  toward  the  newly  heated  region.  On  aver¬ 
age  the  old  channel  moves  toward  the  new  channel.  This  average  motion  does  not  imply  mixing  but 
the  spatial  behavior  of  the  flow  from  the  vortex  filament  pair,  as  will  be  seen  in  the  simulations  of  Sec¬ 
tion  V  (Fig.  16).  effectively  bisects  the  composite  hot  channel  with  a  cold  jet.  The  result  is  two  smaller 
channels  above  and  below  the  original  symmetry  plane  in  which  the  vortex  filaments  are  now  centered. 
A  third  pulse  located  at  v  —  0  between  the  two  modified  channels  will  cause  each  of  these  smaller 
channels  to  bifurcate  again  with  fluid  jets  from  the  top  of  the  upper  channels  and  the  bottom  of  the 
channel  impinging  on  the  third  expanding  channel  from  above  and  below. 


B.  Two-Dimensional  Distortion  (Elliptical  Channel) 


The  above  results  are  readily  generalizable  to  the  other  symmetry  classes,  and  we,  therefore,  will 
not  evaluate  the  form  factors  in  detail  as  we  have  done  above.  To  illustrate  the  effect  of  smooth  distor¬ 
tions  of  a  pulse  (envelope)  from  a  circular  cross  section,  we  will  treat  the  expansion  of  an  elliptical 
pulse.  By  smooth  distortions,  we  mean  in  this  example  that  all  of  the  pressure  and  density  contours  of 
the  pulse  will  be  elliptical  in  shape.  For  a  Bennett  profile,  we  may  express  this  as 


Pix.  y) 


+ 


Pa- P. . _ 

[(.V  -  .v0)2/ft  +  (3<y  -  y0):] 
*<? 


(40) 


In  Eq.  '40),  /?n:  -  aebe ,  where  ae  and  be  are  the  semimajor  and  semiminor  axes  of  the  elliptical 
envelope  of  the  pulse,  and  /3  -  ael be. 


Initially  the  pressure  and  density  gradients  will  be  approximately  parallel.  As  the  channel 
expands,  the  density  distribution  will  retain  an  elliptical  shape  while  the  pressure  distribution  and  the 
flow  pattern  will  approach  radial  symmetry.  Consequently  the  source  term  in  Eq.  (1)  will  be  nonzero 
when  integrated  over  the  expansion  period  r.  In  contrast  to  the  .off  center  beam  propagation  case,  the 
largest  contribution  to  the  vorticitv  comes  from  the  flows  in  the  interior  of  the  expanding  channel 

boundary  <r  ^  R(t)).  since  the  density  is  relatively  uniform  outside.  Our  derivation  will  assume  the 
configuration  of  Fig.  2,  where  the  circle  represents  the  pressure  distribution  and  the  ellipse  forms  the 
envelope  of  the  density  distribution.  Since  we  will  not  account  for  the  noncylindrical  shape  of  the  pres¬ 
sure  distribution  at  early  times,  Eq.  (13)  provides  an  adequate  description  of  the  flow  field  for  this 
model.  The  problem  is  two-dimensional,  as  in  the  case  of  the  first  asymmetry  class;  so  only  is 
nonzero.  The  divergence  of  the  flow  velocity  is  nonzero,  i.e. . 


V  •  v 


2 U 

R 


(41) 
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Because  the  vorticity  generating  flows  are.  therefore ,  compressible,  the  integration  of  Eq.  ( 1 )  is  simplest 
if  we  use  Lagrangian  variables  from  the  beginning.  The  density  of  a  fluid  element  will  vary  it.'  ersely 
with  its  volume  over  time,  as  determined  by  Eq.  (41).  We  may  express  the  average  reduction  ;r  den¬ 
sity  as  the  channel  expands  by  the  equation8 

p(r,  /)  =  pfr.  ■  42) 

where  ^  now  varies  with  time.  R  <r)  is  the  boundary  of  the  outward  flow  as  in  Eq.  ( 13)  and  R0  -  R  <0>. 
Using  Eq.  (42)  and  Eq.  (15)  with  r  -  0,  we  find  that,  to  a  very  good  approximation, 

-rr  In  p(r.  r>  =  In  p(r.  0)  -  -  In  <p„/po>  g(r,  0).  (43) 

d0  ~  dO  ~  oO  ~ 

As  in  Section  IV. A,  we  have  suppressed  the  dependence  on  S0  and.  in  this  case  T0,  to  shorten  our 
notation.  Thus  the  density  distribution  at  very  early  times  dominates  the  density  dependence  of  the 
source  term.  With  Eqs.  (14  and  41-43).  Eq.  (1)  becomes 


2U,  .  P-  d  .  n.  J  U  U2 

nr  *  -  - ln  77  7,  *  ’■■f 


We  may  now  integrate  Eq.  (44)  over  the  time  interval  (0.  r).  As  in  the  case  of  the  density,  the  vorti¬ 
city  present  in  a  fluid  element  at  time  r  <  r  decreases  with  subsequent  expansion  by  the  factor 
R:U)/R:It ).  This  effect  is  associated  with  the  nonzero  divergence  of  the  velocity  and  thus  accounts 
for  the  second  term  on  the  left  hand  side  of  Eq.  (44). 


Our  Lagrangian  variables  are  r  (t )  and  Oir)  -  0  (constant  in  time),  whicn  give  the  position  of  a 
fluid  element  at  ,ime  t  »  r.  From  Eq.  <  1 3 )  we  find  that  Hr)  and  HO)  =  r0  are  related  by 


When  we  express  the  initial  density  distribution  in  terms  of  r(r),  the  vorticity  at  time  r  is  the  integral 
of  the  source  term  in  Eq.  (44)  over  the  time  interval  (0.  r).  reduced  by  the  factor  /? : ( z )//?-< r >  to 
account  for  the  nonzero  divergence  of  the  velocity.  This  gives  us 


(Ar.  0,  r) 


=  -  In 


9  p 

r  R2<r) 

_d_ 

V 

+  iii 

Po 

dO8  R<r)  • 

•'o  R:ir) 

Jt 

R 

R: 

1 _ ,  d  r(r)R  o  „ 

R:U)  °  Po  90 V  R  O)  ■ 


—  ( UR ) 

lit 


(46) 
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Following  the  development  of  Section  A,  we  note  that  the  term  d(UR  )/dt  integrates  to  zero,  and 
we  approximate  the  remaining  portion  of  the  integral  in  terms  of  r  -  tm,  when  the  expansion  flux 
UU)  R';  >  reaches  a  maximum.  This  gives  us 


fAr,  6,t) 


U”T  in 

Poo 

d 

r(r)  R0 

J  n 

2  R2(t) 

Po 

d9  g 

R(t )  ' 

147) 


Comparing  ih.s  with  the  results  of  Section  A,  we  see  that  the  form  of  k.  will  be  quite  similar  in  the  two 
cases. 


We  may  now  integrate  Eq.  (47)  over  the  quarter  plane,  which  contains  a  single  vortex  filament,  as 
shown  in  Fig.  2.  The  integral  is 


-  _  Vi*  , 


Poo 

r“d*  r  ^  8 

z!i. » 

Po 

J  0  Jo  dp 

R 

(48) 


where  r  and  R  are  functions  of  r.  We  may  perform  the  9  integral  and  use  Eq.  (45)  to  transform  from 
r  (t)  to  t)  =  ro/Ro  as  an  integration  variable.  If  we  now  include  the  explicit  dependence  on  and  S0 
we  obtain 


K).  =  — —  In 
2 


Po 


fiAa,  b,  c) 


(49) 


=  Um<R\-  R0)  In 


Pn 


fi-  (a,  b.  c), 


where  a.  b,  and  c  are  given  by  Eq.  (27).  The  form  factor  is 
fiAa.  b.  c)  -  f  dy  y 


giaby,  — .  c)  —  g(aby,  0,  c ) 


(50) 


For  the  initial  effective  radius  of  the  flow  field,  we  use  /?o  ”  aebe.  in  which  ae  and  be  are  the  semima¬ 
jor  and  semiminor  axes  of  the  pressure  pulse. 

We  may  calculate  the  form  factor  f2.  for  a  super-Gaussian  density  distribution. 


—  P  — i 

S0’  ’S0 


exp 


li 

So 


1  + 


cos  29 


(51) 


t 


Inserting  Eq.  (51)  into  Eq.  (50).  we  have 


2-r0  S0 


(52) 


For  S0  =  R o  and  *0  =  1/2  50,  I  f2.\  =  1,  and  our  form  factor  is  a  factor  of  —2  larger  than  for  asym¬ 
metry  class  one.  However,  since  we  did  not  account  for  the  approximate  alignment  of  the  pressure  and 
density  at  early  times,  we  might  expect  ! /2. 1 .  as  given  by  Eq.  (50).  to  be  too  large  by  a  similar  factor. 
Thus,  we  would  expect  a  more  exact  calculation  of  I/2-I  to  yield  values  similar  to  those  of  |/l; I. 


C.  Three-Dimensional  Distortions 


A  pulse  may  also  have  an  axis  which  is  curved  or  kinked.  This  curvature  may  be  the  result  of  a 
perturbation  in  the  envelope  of  the  pulse,  as  in  the  case  of  a  deformable  solid,  or  alternatively,  could 
occur  if  the  deposition  of  energy  along  the  axis  is  nonuniform.  Our  current  treatment  assumes  the 
former  case,  in  which  the  pulse  is  deformed  smoothly,  and  we  will  consequently  parallel  Section  IV. B  in 
our  general  approach.  Figure  3  shows  that  we  choose  the  r-axis  to  coincide  with  the  center  of  mass  of 
the  pulse.  This  causes  the  density  to  vary  as  a  function  of  r  and  0,  and  both  f.  and  are  nonneligible. 
We  will  discuss  each  component  separately,  taking  £.  first. 

C.l  Derivation  of  K). 


The  assumptions  of  Section  IV. B  apply  directly  to  the  derivation  oficj..  since  the  interior  flows  are 
again  responsible  for  the  voriicity  generation.  Figure  3  indicates  that  the  flow  is  radial  from  the  r  axis, 
and  R,y.  the  radius  of  the  expansion  wave,  is  ^  l*0l  S0,  where  *0  ,s  the  displacement  of  the  channel 
axis  from  the  r-axis.  By  convention  ,V0  will  be  positive  for  a  displacement  of  the  channel  to  the  right  of 
the  r-axis.  The  situation  resembles  the  first  asymmetry  class,  off-center  beam  propagation,  although  *0 
oscillates  between  positive  and  negative  values  as  r  changes.  As  in  Section  IV. B.  we  initially  suppress 
r.  and  So  in  our  argument  list  for  .?< s.  r ).  Carrying  over  Eqs.  (41  )-(47).  we  find  that 
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(53) 


in  which  r  and  R  are  functions  of  r.  Because  of  the  similarity  to  asymmetry  class  one.  we  see  that  the 
radial  variable  s  is  given  by  Eq.  (20)  with  ,V0  a  function  r.  and  the  0  integral  is  given  by  Eq.  (21).  We 
may  now  transform  our  radial  integration  variable  from  r(r)  to  v  -  rj R0.  as  was  done  in  solving  Eq. 
(48)  for  asymmetry  class  two.  This  gives  us 


x-Az.r)  = 


Uh 


In 


Poo 

Po 


/j.  6,  t.  :) 


(54) 


where  the  dependence  on  .V0(r  >  and  50  is  now  explicitly  included  through  the  ratios  a.  b.  and  c.  which 
are  defined  in  Eq.  (27).  With  the  assumption  that  g[s,  S0)  ”  Jf  <s/S0),  the  form  factor  is 


p  c  at) 

f)Aa,  b,  c,  r)  *  J  dr\  rj  giabr)  +  c.  z)  -  dr)  ri  gic  -  abr),  r) 

—  f  dr)  r)  g(abr\  —  c,  :). 

»  l  ab 


(55) 


As  an  example  let  us  use  a  Gaussian  function  for  gts/S0),  Eq.  (31).  where  s  is  given  by  Eq.  (20) 

with 


X0(z)  —  h  cos  (2 irz/\) 


(56) 


in  the  coordinate  system  defined  by  Fig.  3.  Equation  (55)  gives  us 
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*o<r)S0 
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(57) 


Note  the  similarity  with  our  calculation  for  /,.  (Eq.  (52)).  We  also  pointed  out  earlier  in  this  subsec¬ 
tion  and  in  Fig.  3  that,  in  our  model  for  the  flow  field,  R0  <  |.Y0i  +  S0,  and,  therefore,  suggest  using 
h  +  S0  as  an  estimate  for  R0.  If  h  =  50,  then  the  average  magnitude  of  X0  is  =2/3  S0  and  we  find 
that  the  average  |/3.|  is  =0.3.  As  for  the  elliptical  channel  case,  the  pressure  will  not  be  symmetric 
about  the  c-axis  at  early  times,  and  our  values  for  |/}.  I  may  be  somewhat  high.  Note  that,  by  Eq.  (56). 
most  of  the  vorticity  is  generated  at  the  ends  and  center  of  the  channel  section. 


C.2.  Derivation  of  Ky, 


To  solve  for  the  0  component  of  the  vortex  filament  strength,  we  return  to  Eq.  G).  Following 
the  line  of  reasoning  which  lead  to  Eq.  (44),  we  find  that 
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Again  the  density  distribution  at  early  times  determines  the  vorticity  through  the  derivative  of 
£<_s.  S0,  0)  in  Eq.  (58).  Proceeding  as  in  subsections  IV. B  and  C.l.  we  use  lagrangian  variables  r(r). 
0,  and  r.  and  integrate  over  time  first  to  obtain 
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As  in  the  case  of  Eq.  (22),  we  simplify  our  integral  by  replacing  H/)  with  the  value  of  r  at  a  specific 
time;  here  we  use  Hr),  which  is  the  most  convenient  choice  for  the  spatial  integral  below.  Again  we 
note  that  d(UR)/dt  integrates  to  zero,  and  we  approximate  the  remaining  portion  of  the  integral  in 
terms  of  t  -  rm.  Thus  we  obtain 
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We  now  integrate  Eq.  (60)  over  the  region  fr«(0,  °°).  r«(0,  A/2)|,  which  contains  vorticitv  of 
one  sign,  in  analogy  to  the  two-dimensional  cases.  The  integral  is 
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We  now  perform  the  r  integral  and  use  Eq.  (45)  to  transform  Hr)  to  rj  =  r(y//?0  as  an  integration  vari¬ 
able.  Including  the  explicit  dependence  on  Xn  and  S0  through  the  ratios  r,y'So  and  Xq/S0  as  before,  we 
have 
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(62) 


where  a.  b,  and  c  are  defined  by  Eq.  (27)  and  the  form  factor  is 

/}„  (a.  b,  c.  8)  “  dr)  rj  8,  A/2,  c)  —  glabr),  8,  0,  c)j.  (63) 

For  the  initial  radius  Rq,  we  may  again  use  LVq!  +  50  If.e<s/S,))  is  given  by  Eqs.  (20).  (31).  and  (56). 
we  find  that  the  form  factor  is 
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Equation  (64)  shows  that  Ik^I  is  maximum  on  the  .v-axis  (8  =  0.  rr).  This  is  analogous  to  the  situa¬ 
tion  for  /3.  which  is  determined  primarily  by  r  —  0,  A/2.  Notice  also  that  the  maximum  values  of  i 
and  \fi»\  agree  at  z  -  0.  A/ 2  and  8  =  0.  rr,  respectively,  as  they  should  because  vorticitv  is  divergence 
free. 


D.  Nonuniform  Deposition  of  Energy 

In  asymmetry  classes  two  and  three  (sections  IV. B  and  C).  we  assumed  that  the  pulse  was  dis¬ 
torted  smoothly,  similar  to  changes  in  shape  of  a  deformable  solid.  Vorticity  generation  also  takes  place 
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when  energy  deposition  is  nonuniform  within  the  -,n "elope  of  a  pulse.  The  simplest  model  of  this  ! 

situation  wou'd  be  that  of  .V  simultaneous  small  pulses  within  the  enveloe  of  the  actual  pulse  and  with 
the  same  integrated  (total)  energy  as  the  actual  pulse.  To  estimate  the  amount  of  vorticity  generated  in 
this  model,  we  may  extrapolate  the  results  of  Section  Iv>  A  for  a  pulse  which  is  misaligned  relative  to  a 

preformed,  hot  channel.  Each  of  the  .V  pulses  will  produce  a  hot  channel  through  which  the  expanding  ] 

shock  waves  produced  by  the  other  .V  -  1  pulses  will  yass.  Taken  separately,  each  pair  of  pulses.  ] 

labeled  by  k  and  /,  would  produce  two  pairs  of  vortex  filaments;  Eq.  (29)  gives  the  approximate 
strength  of  each  filament  kkl.  Given  that  .V  is  a  small  number  and  the  pulses  are  sufficiently  close,  we 
remember  from  Section  IV. A  that  the  strength  kkl  will  vary  slowly  with  the  displacement  of  pulses  k 
and  /.  Thus  we  have  IkwI  =  k  (constant).  Since  there  are  .Y(.V  -  l)/2  pairs  of  pulses,  we  find  that 
convective  mixing  of  characteristic  strength  :V(.V  -  1)  i<  occurs  when  a  small  number  .V  of  simultane¬ 
ous  pulses  are  deposited  close  together.  When  the  displacements  of  the  pulses  are  greater  than  the 

radius  of  the  channel  which  one  of  the  pulses  would  produce  separately  (i.e..  Xtf  >  R llk  =  R\/),  the  A 

form  factor  in  Eq.  (29)  decreases  rapidly  with  increasing  Xq.  Thus  for  a  large  number  of  pulses  or  for 
pulses  which  do  not  overlap,  only  "nearest  neighbors"  will  interact  to  foduce  significant  vorticity.  The 

number  .V*  of  nearest  neighbors  for  each  pulse  will  depend  on  geometry,  and  the  resultant  mixing  ; 

strength  will  be  proportional  to  .V  S'k.  We  expect  that  the  pulse  sizes  will  determine  the  minimum 

scale  of  the  turbulent  structure  observed  in  the  composite  channel  produced  by  the  collection  of  pulses.  | 

We  will  present  further  discussion  of  this  higher  order  application  of  the  theory  and  a  comparison  with 
relevant  experimental  data  in  a  future  article. 

V.  DETAILED  NUMERICAL  SIMULATIONS 

We  have  performed  preliminary  detailed  numerical  simulations  of  this  problem  using  the 
FAST2D9 11  computer  code  in  order  to  validate  and  calibrate  the  approximate  analytic  model  developed 
above.  The  FAST2D  code  solves  the  two-dimensional  equations  for  convervation  of  mass,  momen¬ 
tum.  and  energy,  employing  the  techniques  of  fiux-corrected  transport  <FCT)  and  time  step  splitting 
In  addition  to  accounting  for  shocks  properly,  which  the  theoretical  model  does  not  do.  the  simulations 
are  capable  of  describing  the  late-time  motions  and  profiles  as  modified  by  the  induced  vorticity  Thus 
the  various  elements  of  the  mixing-time  approximations,  and  their  significance,  can  be  evaluated  using 
detailed  calculations  of  representative  problems.  In  this  section  we  will  discuss  simulations  of  examples 
of  the  two-dimensional  asymmetry  classes  (1  and  2).  The  case  of  off  center  pulse  propagation  will 
receive  more  attention,  since  the  vorticity  strength  will  usually  be  greater  than  for  smooth  two- 
dimensional  distortions  and  because  the  general  features  of  the  two  cases  are  the  same  We  have 
addressed  the  question  of  channel  c-  sling  elsewhere'  6  in  the  context  of  convective  cooling  of  lightning 

channels,  and  we  will  discuss  complex  multipulse  turbulent  flows  in  a  future  paper.  I  , 


22 


A.  Off  Center  Pulse  Propagation 


Our  simulations  of  off  center  pulse  propagation  utilize  Cartesian  coordinates  and  the  configuration 
of  Fig.  8.  which  shows  a  stationary,  rectangular  grid  Only  the  upper  half  plane  appears  in  the  calcula¬ 
tion,  since  the  .v-axis  forms  a  symmetr>  axis  The  lower  boundary  is  reflecting  while  the  others  are 
open.  The  grid  consists  of  100  by  30  cells  of  dimension  fix  and  by.  each  varying  from  1  mm  to  5  mm. 
To  resolve  the  channel  dynamics  sufficiently,  we  have  embedded  a  finely  zoned  region  of  50  x  20  cells, 
each  1  mm  on  a  side,  in  the  center  of  the  grid,  where  the  axis  of  an  initial  1.0  cm  radius  hot  channel 
was  positioned.  With  increasing  distance  from  the  central  fine  grid,  the  values  of  5.x  and  by  transitioned 
smoothly  from  1  mm  to  5  mm.  Near  the  sides  and  top  of  the  grid,  5.x  and  by  were  both  5  mm. 

Figure  9  shows  a  plot  of  pi.x.  0.  0)  and  P t.v.  0.  0)  (along  the  y  —  0  plane)  at  the  time  the  second 
pulse  is  deposited.  Bennett  profiles.  1,  1 1  *  s-  Ss,  ):.  for  the  density  and  pressure  deviation  from 
ambient  were  used  in  these  calculations  The  dots  m  on  the  two  curves  show  the  location  and  spacing 
U  mm)  of  the  finite  difference  grid  points  in  the  vicinity  of  the  channel  center. 

To  compute  the  vortex  filament  strength  approximately  we  integrate  the  vorticity  over  the  upper 
half  plane  to  obtain 

x1;(r)  -  J  ^  dx  dy£Ax,  >.  i )  -  J  ^  dxvx(x.  0.  r)  =  ft  ™'dxvx(x.  0.  rl  (65) 

where  ,rmn  and  are  the  respective  limits  of  the  coordinates  on  our  grid.  To  derive  Eq.  (65)  we 
have  assumed  that  v  x ,  v,  —  0  as  x,  y  —  » . 

Figure  10  shows  the  results  of  our  calculations  of  if \.(t)  for  px/po  “  1G,  R0  “  0.4,  R]  —  1.3  cm. 
5 q  *  1.0.  and  .V0  —  1.0  cm  using  ne  density  profile  of  the  channel  formed  by  the  first  pulse  and  the 
pressure  profile  of  the  second  pulse  as  shown  in  Fig,  9.  Throughout  the  simulations,  y  =  cP,  c,  ” 
1.35.  where  cP  and  o  are  the  specific  heats  for  air  with  constant. pressure  and  constant  volume,  respec¬ 
tively.  After  an  initial  positive  transient  and  a  weak  relaxation  oscillation,  the  integrated  vorticity  set¬ 
tles  down  to  a  steady  state  (negative)  value.  We  may  compare  our  formula,  Eq.  (29),  for  the  vortex 
filament  strength  with  the  residual  value  shown  in  Fig.  10.  If  we  assume  Rm  =  1.0  cm.  a  =  0.4.  b  ~ 

1.0.  and  c  -  1.0.  then  /;.  (.4.  1.0.  1.0) - 0.23  from  Fig.  15.  Substituting  these  numbers  into  Eq. 

(29)  with  l’„  =  c,  taken  as  4  x  104  cm/sec  gives  it|.  (r)  =  — 1.9  x  104  cnr/sec.  This  number  is  within 
10%  of  the  value  of -2.1  x  104  cnr/sec  which  we  obtain  from  our  simulations.  The  speed  of  sound  cs 
should  be  evaluated  using  only  the  density  and  temperature  profiles  of  the  preformed  channel  at  the 
location  of  the  center  of  the  second  pulse,  i.e..  x  *  -X0  —  -1  cm  for  the  example  here.  The  charac¬ 
teristic  time  r  must  be  approximately  equal  to  (R ,  -  Rq)/(cs/2)  =  45  n s  to  justify  using  c,  as  the 
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at***. 


effective  maximum  expansion  velocity  of  the  heated  cylinder.  This  value  is  within  25%  of  the  observed 
value  of  r  <  60  >is. 

Using  Eq.  (6)  with  v  -  0  and  our  theoretical  value  for  ic,.  gives  the  maximum  flow  velocity 
expected  in  the  symmetry  plane. 
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(66) 


when  y  —  0  8  cm  (from  the  simulation).  This  value  is  in  good  agreement  with  the  jet  velocity  of  -85 
m/s  measured  in  the  detailed  simulations.  The  mixing  time  calculated  from  Eq.  <10)  using  parameter 
values  from  the  simulation  is  about  360  ps. 


d  v 

In  deriving  Eq.  < 25 )  from  Eq.  <22)  the  time  derivative  term  — ,  which  must  integrate  to  zero 

dt 

d  v 

during  the  expansion,  is  neglected  in  favor  of  the  v~—  term,  which  is  alwavs  of  the  same  sign  and 

dr 

hence  contributes  most  to  the  integral  over  the  interval  (0.  r).  A  posteriori  justification  for  this 
approximation  is  obtained  from  Fig.  10.  which  shows  both  the  general  cancellation  of  a  big  transient 


component,  as  assumed  in  Section  IV.  and  good  quantitative  agreement  with  the  residual  vorticity  term. 
This  means  that  the  amount  of  vorticity  generated  depends  strongly  on  the  eventual  radial  displacement 
of  each  fluid  .lement  but  not  on  the  detailed  expansion  history.  The  scale  of  the  vorticity  is  sei  by  an 
average  expansion  velocity  only 


Figure  1 1  shows  a  schematic  diagram  to  explain  physically  how  the  initial  positive  vorticity  spike 
of  Fig.  10  occurs.  The  expanding  shock,  shown  at  four  different  times  as  a  dashed  line,  travels  faster 
and  drives  the  post  shock  flow  faster  in  the  low  density  material.  We  represent  the  density  gradient  as 
being  localized  between  the  two  solid  semi-circles  in  an  annulus  centered  at  S0.  The  shear  resulting 
from  the  differential  acceleration  of  the  inner  and  the  outer  fluid  corresponds  to  vorticity  of  positive 
sign  w  hich  vanishes,  as  assumed,  when  differential  velocity  ceases.  The  angle  between  Vp  and  V  P 
•  assumed  normal  to  the  dashed  line  shock  fronts)  appears  at  three  locations  along  the  original  channel 
periphery.  The  residual  vorticity.  which  corresponds  to  that  appearing  in  Eqs.  <25)  and  < 28 > .  in  effect 
measures  the  net  outward  expansion  of  the  new  channel. 

Figures  12  and  13  show  the  vorticity  strength  ic,_. <r>  for  several  values  of  .V,i  with  the  other 
remaining  variables  in  our  standard  case  held  constant.  When  the  axis  of  the  second  pulse  fell  within 
5,)  of  the  center  of  the  channel  •  Fig.  12).  the  transient  term  appeared  to  oscillate  more  rapidly  than  the 
cases  for  which  the  axis  of  the  second  pulse  fell  "outside"  the  original  channel  (Fig.  13)  The  residual 
values  it |..  were  also  more  nearly  equal  for  the  cases  in  Fig.  12  than  for  those  in  Fig.  13.  The  general 
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trend  of  <<,.(/  —  °o)  in  the  five  cases  of  Figs.  12  and  13  follows  qualitatively  the  shape  o  f  f,Ja.  b,  c) 
plotted  for  various  cases  in  Figs.  5,  6,  and  7.  The  uncanceiled  vorticity,  which  we  called  kI:  in  Section 
Ill.  decreases  from  a  maximum  around  .V0  =  S()  both  as  .V,)  —  0  and  as  X,)  becomes  large  compared  to 
The  ratio  of  kl.  at  .V0  ”  0.2  cm  to  kl.  at  ,T0  “  0.1  cm  is  probably  not  as  small  as  expected  from 
the  curves  J\.  in  Fig.  6  because  the  sound  speed  is  higher  when  the  second  pulse  is  nearer  the  center  of 
the  preformed  channel.  The  increasing  sound  speed  partially  counterbalances  the  tendency  of  /,.  to  fall 
off  as  .V0  approaches  zero. 


The  major  issue  in  interpreting  these  calculations  using  the  formulae  of  Section  III  revolves 
around  the  choice  of  U^t/2  for  Eq.  (28).  For  estimates  in  the  case  of  off  center  pulse  propagation,  w  - 
have  assumed  L„  is  the  sound  speed  of  the  fluid  at  the  center  of  the  second  pulse  just  before  the 
second  pulse  is  actually  deposited.  We  have  also  assumed  that  the  average  expansion  velocity  is  LiJ 2. 


Figure  14  compares  the  effects  of  initial  pulse  size  on  the  second  pulse,  given  constant  values  of 
the  channel  parameters  and  the  overpressure  and  displacement,  ,Vrt,  of  the  second  pulse.  We  notice 
that  the  transient  term  appears  to  scale  in  magnitude  as  Ra  while  the  value  of  the  "residual"  vorticity 
appears  to  scale  as  /?02-  We  may  explain  the  scaling  of  the  residual  vorticity  in  terms  of  our  analytic 
"late  time"  treatment.  Instead  of  using  the  estimate 

Vm  =  Cj  (67) 


we  may  calculate  L’m  from  the  equation 


U 


20? |  -  R„) 

T 


(68) 


where  r  is  the  time  of  expansion  from  R0  ^i-  1°  P'S-  14,  we  see  that  the  extrema  in  for  the 
cases  of  R0  —  0.2  cm  and  0.4  cm  occur  at  almost  the  same  values  of  t  and  that  the  relaxation  to  the 
final  value  of  *1;  occurs  at  the  same  rate  for  the  two  cases.  Thus.we  have  (for  any  consistent  definition 
of  r)  the  relation 
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(69) 


However  we  expect  that 


! 


(70) 
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Thus  Eq.  (68)  yields 


Um<  =  j  UmB.  (72) 

Estimates  of  the  form  factor  indicate  that 

ft:  =  ft. 

Using  Eqs.  (31)  and  (46) - (5 1 ) ,  we  find  that 

*  I;  __  ~  RqA)‘  __  _1_ 

k*  '  <*18-  Rob):  "  4 

for  the  case  shown  in  Fig.  14. 

This  analysis  indicates  that  c,  provides  only  a  crude  estimate  of  L'm  as  the  energy  deposited 
becomes  small.  From  Eq.  (69)  and  the  fact  that  only  the  preformed  channel  was  identical  in  cases  -1 
and  8,  we  conclude  that  the  relaxation  time  of  the  transients  depends  as  much  on  the  already  existing 
density  gradients  as  it  does  on  the  expansion  of  the  new  channel.  The  ambient  sound  speed  at  ,V0  - 
1.0  cm  is  about  4  x  104  cm/sec  for  the  standard  case  and  the  reduced  energy  case  of  Fig.  14.  The  final 
relaxation  to  the  long-time  value  x1;  begins  at  about  40  ps  for  both  cases,  about  the  time  required  for  a 
sound  wave  to  reach  the  far  side  of  the  original  hot  channel  from  the  edge  of  the  second  pulse  (a  dis¬ 
tance  of  1.8  cm  and  1.6  cm  for  cases  A  and  B  respectively). 

In  Table  1  and  Fig.  15,  we  present  estimates  of  the  form  factor  /,.  for  the  case  S0  -  1.0  cm.  Rn 
"  0.4  cm,  /?|  —  1.4  cm,  p»/po  “  10  (channel  formed  by  first  pulse)  and  PJ Px  -  31  (second  pulse), 
as  shown  in  Figs.  8  and  9.  These  values  correspond  to  a  =  0.4  and  b  ==  1.0.  We  vary  c  from  0.2  to 
2.0.  To  estimate  the  value  of  Um,  we  used  Eqs.  (67)  and  (68).  For  Eq.  (68).  we  define  r  to  be  the 
time  elapsed  between  the  deposition  of  the  second  pulse  (t  -  0  in  Figs.  13  and  14)  and  the  last 
extremum  in  St.  before  the  monotonic  relaxation  to  the  residual  value.  The  solid  curve  is  the  estimate 
from  Eq.  (30)  using  a  Bennett  density  profile,  Eq.  (32),  for  the  original  channel.  Near  c  —  1.0.  Eq. 
(67)  appears  to  provide  a  better  estimate  of  im  than  does  Eq.  (68).  at  least  when  R 0  —  1.0  cm.  Better 
agreement  between  the  values  of  /,.  derived  from  simulations  and  the  theory  is  expected  when  the 
energy  deposition  is  slow  and  the  adiabatic  expansion  treatment  is  a  better  approximation. 

Finally,  to  illustrate  the  details  of  the  flow  resulting  from  off  center  •  ise  propagation,  we  have 
performed  a  simulation  of  two  sequential,  noncollinear  pulses  without  the  use  of  a  symmetry  plane. 
The  Cartesian  grid  consists  of  100  x  100  cells  with  50  x  50  uniform  fine  cells  in  the  center.  The 
dimensions  of  the  remaining  cells  increase  geometrically  to  move  the  boundaries  far  away  from  the  hot 
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channel.  The  first  pulse  deposits  energy  in  the  center  at  time  r  -  0.0  while  the  second  occurs  1  ms 
later  and  deposits  energy  to  the  right  of  center  ( XQ  -  -0.5  cm).  The  peak  overpressures  (/>„,„  ~  /*») 
at  the  instants  when  energy  deposition  occurs  are  4.7  atm  and  2.4  atm,  respectively.  Again  the  pulses 
have  Bennett  profiles.  Eq.  (32),  with  a  Bennett  radius  S0  *  0.5  cm;  thus  we  have  c  -  1.0.  Figure  16 
shows  density  contours  and  velocity  vectors  at  t  —  1.26  ms  and  density  contours  at  t  —  2  ms  for  the 
finely  gridded  central  region.  The  noncollinearity  of  the  pulses  has  produced  an  elongated,  hot  central 
core.  In  addition,  the  vector  velocity  plot  clearly  shows  the  enhanced  flow  between  the  two  vortex 
filaments,  which  are  separated  by  2v  ~  1.1  cm  or  approximately  2  Bennett  radii.  This  flow  pulls  the 
original  hot  channel  (produced  by  the  first  pulse)  toward  the  center  of  the  second  pulse,  as  predicted  by 
our  analytic  theory.  The  diagram  corresponding  to  t  -  2  ms  shows  that  the  center  of  the  original  chan¬ 
nel  has  cooled  (has  a  higher  density)  as  a  result  of  this  displacement.  We  observed  identical 
phenomena  in  all  similar  simulations. 

B.  Two-Dimensional  Distortions 

To  simulate  vorticity  generation  by  smooth  two-dimensional  distortions,  we  chose  the  case  of  an 
elliptical  pulse  as  defined  in  Fig.  17.  The  parameters  in  Fig.  17  are  typical  of  elliptical  laser  pulses  used 
in  recent  studies  by  Greig  et  al.i:  We  note  that  the  pulse  is  much  weaker  than  those  considered  in  Sub¬ 
section  IV'  A,  and  thus  expect  the  residual  vorticity  to  be  lower  in  value.  We  use  a  constant  value  of 
y  -  1.4  for  the  ratio  of  principal  specific  heats.  The  Cartesian  grid  is  quite  similar  to  that  used  in  the 
previous  case,  although  more  cells  (150  x  75)  are  used.  The  embedded  fine  grid  consists  of  100  x  50 
cells,  and  outside  that  region  the  x  and  >  dimensions  of  the  cells  increase  geometrically.  This  moves 
the  boundaries  far  from  the  region  of  interest,  permitting  the  shock  to  propagate  well  away  from  the 
hot  channel  before  reaching  the  boundary  After  pressure  equilibrium  of  the  hot  elliptical  channel  is 
reached,  we  numerically  elimina  e  the  radial  flows  associated  with  the  shock  and  follow  only  the 
incompressible  flows  associated  with  cor  mechanism  for  channel  cooling.  This  permits  use  of  longer 
time  steps  (factor  of  — 100)  and  reduces  the  running  time  significantly. 

Figure  18  shows  a  sequence  of  density  contour  diagrams  for  this  calculation.  Defining  the  time 
/  —  0  to  coincide  with  the  instantaneous  deposition  of  energy,  we  see  the  shock  wave  moving  out  of 
the  fine  grid  and  leaving  a  hot  channel  with  an  elliptical  cross  section  at  t  ~  50  ms.  The  channel 
remains  approximately  the  same  for  /  =  1  ms  and  begins  to  distort  noticeably  by  t  ~  2  ms  under  the 
influence  of  the  flow  shown  in  Fig.  2  This  flow  pattern  continues  to  pull  the  channel  outward  along 
the  x-axis  as  the  two  vortex  filament  pairs  apparently  move  apart.  Because  a  relatively  coarse  grid  is 
used,  these  calculations  do  not  resolve  smaller  scale  turbulent  structure  which  has  been  observed  exper¬ 
imentally.12  In  addition,  perturbations  introduced  experimentally  might  alter  the  interactions  among 
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vortex  filaments  and,  therefore,  the  flow  pattern  shown  in  our  calculations  at  late  times.  However,  this 
simulation  should  provide  a  useful  treatment  of  the  larger  scale  flows  which  determine  the  rate  at  which 
the  ambient  atmosphere  mixes  with  the  hot  channel  gas  as  a  result  of  smooth  two-dimensional  distor¬ 
tions  of  a  pulse. 

To  obtain  an  estimate  for  the  strength  of  each  vortex  filament,  of  which  four  are  present,  we 
integrate  the  vorticity  over  the  upper  right  quadrant  (xt(0,  =®  >.  v « < 0.  °°)|  to  obtain 

K-.(r)  —  f  dx  v.(x.  0,  r)  —  f  dv  v>  (0.  v.  t)  (75) 

-  J  a  '  Jo  • 

J*  'max  .  ,  C  k  max  .  _ 

dx  v.(x,  0.  t)  -  j  dv  v.  <0.  v,  t ) 

where  xm„  and  vma,  are  the  positive  limits  of  our  grid.  In  deriving  Eq.  (75)  we  have  again  assumed 
that  vt,  v,  —  0  and  x.  y  —  °°.  Figure  19  shows  the  evolution  of  k:.  for  the  elliptical  channel  simula¬ 
tion.  As  we  have  indicated  in  Fig.  2,  the  residual  vorticity  in  the  upper  right  hand  quadrant  should  be 
positive,  the  flow  being  inward  at  the  ends  of  the  major  axis  and  outward  along  the  minor  axis.  During 
the  expansion  to  pressure  equilibrium,  however,  the  acceleration  term  dv/3r  induces  a  transient  nega¬ 
tive  vorticity.  By  our  previous  definition,  rmix.  appears  to  be  <75  ps  and  the  residual  vorticity  is  —2.7 
x  102  cm:/s.  Defining  R:(i)  -  ait)  bit)  where  alt)  and  bit)  are  the  semimajor  and  semiminor  axes 
of  the  ellipse,  we  have  R0  —  0.25  cm.  Adiabatic  expansion  would  give  us  Ri  —  0.36  and  p^/po  ~2 
for  the  hot  channel  at  pressure  equilibrium,  and  from  the  simulation,  we  obtain  Um  <  3  x  103  cm/s. 
which  is  less  than  the  ambient  value  of  c's  by  an  order  of  magnitude  From  Figs.  2  and  17,  we  find  that 
5o  =  0.31  and  X0  —  0.19.  From  Eq.  (52)  and  the  ensuing  discussion,  we  obtain  |/2.|  =  1.  We  may 
now  use  Eq.  (49)  to  find  an  estimate  of  |k2.|  (theory)  <  2.3  x  10:  cm2/s.  Theory  and  simulation, 
therefore,  agree  within  a  factor  of  20%. 

VI.  CONCLUSIONS 

This  paper  represents  a  first  step  in  analyzing  the  cooling  of  hot  channels  produced  by  local  energy 
ueposition  in  a  gaseous  medium.  We  have  identified  a  mechanism  which  produces  sufficient  rotational 
motion  to  account  for  the  cooling  rates  of  electric  discharge  channels  produced  in  the  laboratory.'-5  As 
the  channel  expands  to  pressure  equilibrium,  a  shock  is  produced,  and  asymmetries  between  pressure 
and  density  gradients  generate  vorticity  according  to  Eq.  (1).  After  pressure  equilibration  and  propaga¬ 
tion  of  the  shock  wave  well  away  from  the  hot  channel,  vorticity  is  no  longer  generated;  however, 
significant  residual  vorticity  exists.  We  may  represent  the  resulting  flow  field  in  terms  of  one  or  more 
vortex  filament  pairs  of  strength  ±  »c,a,  where  i  -  1,  2.  3  labels  the  type  of  asymmetry  and  a -corn, 
labels  the  vector  component  of  the  vorticity.  The  strength  *  satisfies  the  equation 
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where  Um  is  the  expansion  velocity  of  the  channel  boundary  when  the  expansion  flux  U(t)  R(t)  is  a 
maximum;  p*  and  p0  are  the  ambient  density  and  the  density  at  the  center  of  the  hot  channel  just  after 
energy  deposition;  and  /,„  is  a  form  factor  which  may  be  evaluated  from  the  formulae  presented  in  Sec¬ 
tion  IV.  For  asymmetry  class  1  (off  center  pulse  propagation),  Um  =  c\  since  the  flows  at  the  channel 
boundary  generate  most  of  the  vorticity.  For  smooth  distortions  of  the  pulse  envelope,  which  belong  to 
asymmetry  classes  2  and  3,  Um  <  cs,  since  the  subsonic  interior  flows  generate  most  of  the  residual 
vorticity.  We  have  also  identified  another  subclass  of  pulse  distortions  which  arise  from  energy  deposi¬ 
tion  (within  the  envelope  of  a  pulse)  which  does  not  follow  a  smooth  functional  form  such  as  the  Ben¬ 
nett  profile  of  Eq.  (40).  As  a  simple  model  we  might  consider  deposition  to  occur  as  several  simultane¬ 
ous  packets  of  energy  within  the  boundary  of  the  pulse.  In  a  future  paper  we  will  analyze  and  simulate 
this  situation  in  detail  using  the  FAST2D  code.  Accurate  numerical  simulations  of  examples  of  the 
two-dimensional  asymmetry  classes  have  revealed  good  agreement  with  Eq.  (76). 


With  an  estimate  for  y,  which  appears  from  simulations  to  be  =  S,)  and  Eq.  (76)  for  x,a.  we  may 
obtain  values  for  the  mixing  time  scales  rmix  from  Eqs.  (lOa-e).  For  the  two-dimensional  asymmetry 
classes  7mix  is  the  approximate  time  for  doubling  of  the  volume.  For  three-dimensional  distortions,  Kit> 
and  Ki;  are  both  nonnegligible,  and  the  cumulative  effects  of  the  two  components  must  be  calculated  to 
obtain  a  reasonable  doubling  time  for  the  channel  volume.  The  values  of  rmix(:)  and  rmjx(s>  given  by 
Eqs.  (10a)  and  (10e)  do,  however,  provide  a  useful  estimate  of  the  importance  of  each  component  rela¬ 
tive  to  the  other  and  to  other  types  of  asymmetry. 
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Table  1  —  Values  of  from  numerical  simulations  The  displacement  of  the  axes  of  the  hot  channel 
and  the  latest  pulse  ir  =  V,  v„>  varies  from  0. 2-2.0.  The  table  also  compares  two  methods  of  computing 
The  sound  speed  is  computed  at  the  position  of  the  pulse  just  prior  to  energy  deposition. 
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Fig.  I  —  Schematic  of  asymmetry  class  one.  m  which  a  cylindrical  pulse  of  radius  R ,  propagates  parallel  to  a  hot  channel  of 
h.tracteristic  radius  S«  formed  by  a  previous  pulse.  The  second  pulse,  offset  from  the  first  pulse  by  V,,.  expands  to  radius  R , 
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Ft*  -  —  The  second  asvmmetrv  class  includes  pulses  with  smooth,  two-dimensional  deviations  from  a  circular  cross  section 
This  schematic  shows  the  upper  hall  ol  the  flow  held  produced  by  a  pulse  with  elliptical  energy  contours.  After  deposition,  the 
pressure  distribution  and  velocity  tie  id  approach  cylindrical  symmetry  while  the  density  distribution  remains  elliptical.  We  model 
this  analytically  by  assuming  a  cylindrical  expansion  from  radius  R,  to  R  :n  the  presence  of  an  elliptical  density  distribution 
characterized  by  a  radius  .Vp.  the  asymmetry  parameter  X \>.  a  density  />,,  at  the  center,  and  an  amhent  density  iix  The  curved 
arrows  indicate  the  direction  ot  residual  rotational  flows  corresponding  to  the  two  vortex  filament  pairs  ol  strength  =  pro¬ 
duced  by  the  asymmetry 
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Ftg  3  -  VV-  model  jsvmmeiry  class  three  by  treating  a  tinu.soidjl  channel  section  of  *a\eiengib  \  This  will  generate  nonnegi 
eible  :  and  M  components  of  the  v orticttv  The  tlo*  field  approaches  cylindrical  symmetr>  about  the  center  of  mass  t.-a\i?>i  vvhi 
the  density  JiNtnbuuon  is  offset  from  the  e-j\i>  by  V».  *hich  \anes  ^ith  position  along  the  ;-a\is. 


Tmix 


Fig  4  —  Sxmmemc  distribution  of  localised  xoructtx  appears  as  two  extended  xortex  filaments  of  strength  —  i<  and  separation  Jr 
The  How  xclov.it>  induce.!  hx  these  xortices  is  shown  along  the  .v-uxis  halfwax  between  the  filaments  If  we  huxe  a  hot  channel  of 
radius  about  which  the  iilamcnts  are  centered,  we  max  define  a  mixing  time  scale  as  the  'mie  required  for  a  fluid  element  at 
v  —  -.S  i  to  reach  v  -  S.,  We  exclude  the  effects  of  smaller  scale  turbulent  motion  which  might  exist  at  the  boundaries  or  in  the 
interior  of  the  xbannel  Here  we  acknowledge  the  existence  of  j  similar  figure,  appearing  in  Journal  of  ihe  Ximospheric  Sci¬ 
ences  '  i 
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C  - ► 

hi  5  —  Integrated  'nrticil>  form  factors  assuming  a  -urcr-Ciaussian'  densitt  pr>>riie 
Laraer  variation  occur-)  with  h  in  the  upper  figure  than  with  j  in  the  lower  tigurc  ", 
tind  peak  corticate  generation  in  the  vicimtc  nt  \  *  ft  ,  tor  til  cases  The  maximum 
maenituJe  ■>(  the  form  factor  ranees  from  11  1  to  0  4  with  11 being  a  reprcserttatu e  '.due 
lor  uuick  estimates 
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Fig.  F  —  Imegra 

led  vorticitv  form  labors  assun 

vng  a  Ben 

ineu  JenMiy  nrotV.e 

:o-n  ihi 

form  faeiors  var 

y  more  \*tth  *  than  vuth  j  No 

te  ihai  api 

preciahie  form  factors  ea 

in  e\is 

a!  V.,  =  O  ’  .V,i  « 

*  hen  n  is  small. 

—  Comparison  of  iniegrated  vorticitv  form  factors  for  three  distinct  density  profiles  at  a  —  h  m  i  %  2  The  variation 
among  form  factors  is  relatively  small  for  widely  differing  density  profiles. 
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F'g  S  —  Initial  conditions  lor  numerical  simulations  ol'  \ orticits  generation  when  a  puise  propjgates  parallel  to  hut  offset  from  in 
existing  hot  channel  The  simulations  include  onl>  the  upper  hall'  plane  with  a  reflecting  lower  boundary  and  outflow  conditions 
at  the  other  boundaries.  By  stretching  the  outermost  /ones  by  a  factor  of  live,  we  mine  the  boundaries  lar  Irom  the  active  re- 
gion 


Fig  d  —  Pressure  and  density  profiles  along  the  symmetry  t  v-r>  plane  when  the  axes  of  the  hot  channel  and  the  latest  pulse  arc 
separated  by  I  0  cm  The  characteristic  radius  of  the  channel  is  10  cm  and  that  of  the  latest  pulse  is  0  4  cm  'Figs  I  and  Si  The 
data  points  correspond  to  cell  centers  and  indicate  the  resolution  provided  by  die  line  /ones. 
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S0  =  1.0  cm 
R0  =  0.4  cm 
X0  =  1.0  cm 

Rnox/^o  =  3 1 

Po.  /P0  =10 

*iz  =  -  2.1  x  104  cm2/s 
Vj,t  =  _85  m/s 
y  =  0.8  cm 

"^mlx  =  360^ts 
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Fig  10  —  Time  history  of  the  vorttcuv  strength  in  the  upper  halt  plane  for  our  standard  case  The  integrated  corn¬ 
et  ty  over  the  entire  plane  is  eero.  since  an  oppositely  directed  vortex  lilanient  having  a  strength  m  euu.ii  magnitude 
is  generated  in  the  lower  half  plane 


Fig  II  —  Scherruiic  diagram  showing  how  the  initial  xorticity  transient  arises  physically  tor  the  situation  in  Fig  I  The  euannty 
<b  is  the  angle  between  the  pressure  gradient  of  the  expansion  llow  and  the  density  gradient  ol  ihe  initial  not  cnunnel  When  the 
outwardly  expanding  lluid  eventually  decelerates  to  rest,  the  acceleration  'and  therefore  the  pressure  gradient'  must  change  sign, 
causing  most  of  the  transient  vorticity  to  cancel 
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Fig  14  —  Integrated  vorticuy  k lot  two  cases  differing  only  <n  the  radius  of  the  second  pulse  Curse  B  is  the  standard  case  and 
curve  A  has  one-fourth  the  energy  per  unit  length  in  one-fourth  the  area,  so  that  the  peak  overpressures  are  equal  The  tran¬ 
sients  differ  in  amplitude  h>  a  factor  of  two  but  the  residual  integrated  vorticuy  differs  by  a  factor  of  four 


Fig  16  —  ' j >  Density  contours  and  t b )  velocity  vectors  at  a  time  ;  «•  I  26  ms  and  (cf  density  contours  at  r  ■*  2  00  ms  Tor  two 
sequential  noncnlhnear  pulses.  The  tirst  pulse  occurred  at  r  *■  0  0  ms  in  the  center  of  the  grid,  and  the  second  pulse  occurred  at 
r  =  I  0  ms  with  a  displacement  of  1)  5  cm  to  the  'reader’s!  right  of  center  The  figures  correspond  to  the  tine  grid  The  density 
contours  vary  from  3  <)  x  Id-4  (innermost'  to  I  I  x  |tr'em/cmJ  (outermost!  In  (hi.  we  have  plotted  a  velocity  vector  at 
every  other  cell  center,  thus  the  centers  of  the  two  vnrtev  filaments  are  somewhat  indistinct  The  major  flow  occurs  between  the 
filaments  and  the  direction  of  ihe  flow  is  from  (he  center  of  the  earlier  pulse  to  the  center  of  the  later  pulse. 


Fig  16  —  ( a >  Density  contours  and  <b)  velocity  vectors  at  a  time  r  «  I  26  ms  and  <c)  density  contours  at  /  =*  2.00  ms  lor  two 
seuuential  noncollinear  pulses.  The  first  pulse  occurred  at  r  -  0.0  ms  in  the  center  of  the  grid,  and  the  second  pulse  occurred  at 
r  =  10  ms  with  a  displacement  of  0  5  cm  to  the  <reader's>  right  of  center  The  figures  correspond  to  the  tine  grid  The  density 
contours  vary  from  2  0  x  10— 1  iinnermosl)  to  I  1  x  10“'  govern3  'outermost)  In  (b>.  we  have  plotted  a  velocity  vector  at 
every  other  cell  center,  thus  the  centers  of  the  two  vortex  filaments  are  somewhat  indistinct.  The  maior  flow  occurs  between  the 


filaments  and  the  direction  of  the  flow  is  from  the  center  of  the  earlier  pulse  to  the  center  of  the  later  pulse 
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Fig  lb  —  i j )  Density  contours  end  'h'  velocity  vectors  ji  a  time  r  *  1.2b  ms  and  'c >  density  contours  at  t  **  2  do  ms  tor  two 
seuuentiai  noncoihnear  pulses  The  first  pulse  occurred  at  r  ”  0  0  ms  in  the  center  of  the  grid,  and  the  second  pulse  occurred  a; 
;  =  I  0  ms  with  a  displacement  of  0  5  cm  to  the  'reader's!  right  of  center  The  figures  correspond  to  the  fine  grid  The  density 
contours  vary  from  .VO  x  I  ft-4  'innermost!  to  I  1  x  ICT1  gm/cm-  foutermosti  In  t b > .  we  have  plotted  a  velocity  vector  at 
every  other  cell  center;  thus  the  centers  of  the  two  vortex  filaments  are  somewhat  indistinct.  The  major  flow  occurs  between  the 
filaments  and  the  direction  of  the  flow  is  from  the  center  of  the  earlier  pulse  to  the  center  of  the  later  pulse. 


Fig  I'  —  Tins  schematic  shows  the  initial  conditions  jnd  Cartesian  grid  lor  our  numerical  simulation  of  voruetiy  generation  hi 
an  elliptical  pulse  We  simulate  only  the  portion  of  the  pulse  above  tne  v-r  symmetry  plane,  surrounding  the  pulse  with  a  uni¬ 
form  grid  of  100  x  50  tine  cells.  The  dimensions  of  the  remaining  cells  then  increase  geometrically  with  distance  from  the  Pulse 


Fig.  18  —  Results  of  a  computer  simulation  of  the  time  deielopment  of  a  laser  pulse  haung  smooth  elliptical  distortion 
The  simulation  used  the  initial  conditions  described  in  Fig  1 7  The  density  contour  diagrams  shown  here  correspond  to  the 
following  time  intervals  elapsed  from  the  point  of  instantaneous  energy  deposition  t a  >  5b  yis.  ' b •  "5  us.  <c*  2  0  ms.  < d  •  4.0 
ms.  and  let  6  4  ms  The  sis  equally  spaced  contour  values  range  from  "4  x  |0~J  g/errr  'contour  1'  to  I  2b  x  !0~'  gem' 
'contour  hi.  the  contour  values  increase  with  increasing  vertical  displacement  from  the  center  of  the  channel.  In  (a*,  the 
shock  appears  just  inside  the  tine  grid  Diagram  'hi  shows  contours  1-5  and  the  others  have  only  contours  2-5.  indicating 
that  the  channel  gas  has  cooled  somewhat  The  flattened  tops  of  the  contours  result  from  clipping  by  the  flux-corrected 
transport  algorithm 
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Fig.  18  —  Results  of  .1  computer  simulution  of  the  time  development  of  j  User  pulse  having  smooth  elliptical  distortion 
The  simulation  used  the  initial  conditions  described  in  Fig  1"  The  density  contour  diagrams  shown  here  correspond  to  the 
following  time  intervals  elapsed  from  the  point  of  instantaneous  energy  deposition  (a)  .'6  /as.  'hi  "5  ms.  <cl  2.0  ms.  ' d >  4  0 
ms.  and  ie 1  b  4  ms.  The  six  equally  spaced  contour  values  range  from  "4  *  |0-J  g/cnv‘  'contour  1 1  to  I  2b  *  10-‘  g,  cm-’ 
'contour  bi.  the  contour  values  increase  with  increasing  vertical  displacement  from  the  center  of  the  channel  In  'a-,  the 
shock  appears  lust  inside  the  tine  grid  Diagram  ' b >  shows  contours  1-5  and  the  others  have  only  contours  2-5.  indicating 
that  the  channel  gas  has  cooled  somewhat.  The  flattened  tups  of  the  contours  result  from  clipping  ny  the  llux-corrected 
transport  algorithm 


Fig  18  —  Results  of  a  computer  simulation  of  the  time  development  of  a  laser  nuisc  having  smooth  elliptical  distortion 
The  simulation  used  the  initial  conditions  described  in  Fig  I"  The  density  contour  diagrams  shown  here  correspond  to  the 
following  time  intervals  elapsed  from  the  point  of  instantaneous  energy  deposition  'a'  sb  o'-  'b1  ”>  u'-  'C  2  '1  ms.  ‘ J '  4  0 
ms.  and  '  e  *  6,4  ms.  The  six  equally  spaced  contour  values  range  from  *4  x  l'l-1  g  vm  'contour  I'  to  !  2b  «  Id"  g,  cm  ’ 
'contour  61:  the  contour  values  increase  with  increasing  vertical  displacement  from  in  center  -t  'he  channel  In  'a',  the 
shock  appears  just  inside  the  tine  grid  Diagram  'b'  shows  contours  1-5  and  the  others  have  niv  contours  2-.v  indicating 
that  the  channel  gas  has  cooled  somewhat.  The  flattened  tops  of  the  contours  result  trom  dirring  Ki  the  !luv-correcied 
transport  algorithm 
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Fic  18  —  Results  oi'  a  computer  simulation  of  the  time  development  of  a  laser  pulse  having  smooth  elliptical  distortion 
The  simulation  used  the  initial  conditions  described  in  Fig.  I'  The  densitv  contour  diagrams  shown  here  correspond  to  the 
following  time  intervals  elapsed  from  the  point  of  instantaneous  energv  deposition  t a )  'ftps.  < b <  "5  pS.  <c'  2.0  ms.  < d <  4  0 
ms.  and  i e '  t>  4  ms  The  six  eduallv  spaced  contour  values  range  from  *  4  x  It)"1  g/cm-  (contour  I  >  to  1  2b  x  10“  '  g,  cm- 
'contuur  ft>.  the  contour  values  increase  with  increasing  vertical  displacement  from  the  center  of  the  channel  In  1  a i .  the 
shock  appears  just  inside  the  hne  grid  Diagram  'b*  shows  contours  1-5  and  the  others  have  onlv  contours  2-5.  indicating 
that  the  channel  gas  nas  cooled  somewhat  The  flattened  tops  of  the  contours  result  from  clipping  b>  the  flux -corrected 
transport  algorithm 
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Fig  18  —  Results  of  a  computer  simulation  of  the  lime  development  of  a  laser  pulse  having  smooth  elliptical  distortion 
The  simulation  used  the  initial  conditions  described  in  Fig  1*  The  densitv  contour  diagrams  shown  here  correspond  to  the 
following  time  intervals  elapsed  from  the  point  of  instantaneous  energv  deposition  1  a »  an  ps.  'h>  "5  ps.  ' c '  2  0  ms.  'd>  4  0 
ms.  and  *e>  ft. 4  ms  The  six  equal).'  spaced  contour  values  range  from  '  4  x  |0“'lg,cm  '  (contour  I  •  to  I  ’t>  x  In"  g  em 
(contour  ft>.  the  contour  vaiues  increase  with  increasing  vertical  displacement  from  the  center  of  the  channel  In  <a • .  the 
shock  appears  just  inside  the  fine  grid  Diagram  < b >  shows  <.onlourt  and  the  others  have  nnh  contours  2-5  indicating 
that  the  channel  gas  has  cooled  somewhat  The  flattened  tops  of  the  contours  result  from  clipping  hv  the  flux-corrected 
tr  msport  algorithm 
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